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ABSTRACT 


A model of the error detection process for the 
testing of software has been developed to investigate 
the relationship between computer program structure 
and error detection and test effort. The model has 
been implemented as a simulation. 
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I. INTRODUCTION 


Software costs are an increasing fraction of total 
computer project costs and, for many computer projects, 
software costs dominate hardware costs. Much of the 
software development costs are for testing, debugging and 
integration; a significant part of the costs after releasing 
the software are for correcting errors. Thus there is 
current interest in the error characteristics; number, type 
(overflow, seguence control) and location of software errors 
in a program. It is generally accepted that computer 
programs with a complex structure* that is one with a high 
incidence of branch instructions and loops, are harder to 
debug and test and more errors persist after release than 
for programs with a more simple structure. An error 
simulation model 1 is presented here which investigates the 
relationship of program structure to error detection and 
test effort. 


Since structure can be controlled during the design 
phase and measured through all phases of a computer project, 
the study of the relationship between structure and error 
characteristics is valuable to the manager of a software 
project. Complex program structures with poor error 
characteristics should be avoided. Poor error 
characteristics result, when many errors are located in 


A The suggestion to 
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consultant to the Naval 
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complex structures in such a way that error detection would 
prove difficult during testing. In cases where complex 
program structures may be necessary to help meet program 
size or speed limitations, it is useful to have an 
indication of the additional testing which may be caused by 
complex structures. it is also useful to be able to compare 
the error characteristics of design alternatives that have 
different program structures. 




* 
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II. NEED FOR RELIABLE SOFTWARE 


A. SOFTWARE COSTS 


The cost of computer systems can be divided into 
hardware (cost of computers and peripheral equipment) and 
software (cost to design, write, test and maintain computer 
programs). The proportion of total costs that is 
attributable to software has been increasing, projections of 
higher costs for people and lower costs for hardware 
components being a continuing trend. The direct cost of 
software is enormous; one estimate of the cost in the United 
States, according to Ref. 1, is 10 billion dollars per year. 
Indirect costs due to late delivery, nonperformance due to 
errors, wrong actions due to erroneous output, etc., are 
also very large. 

Although there is not extensive data available on the 
cost of software projects, there is enough data to indicate 
the magnitude of software costs. References 1 and 2 contain 
a comprehensive study of Air Force computer costs. The Air 
Force estimates that software costs in 1972 were from 1 to 
1.5 billion dollars, which is 4 to 5 percent of the total 
Air Force budget, while computer hardware costs were from 
300 to 400 million dollars. The percentage of computer 
costs due to software has increased from 30 percent in the 
late 1950's to 70 percent now and the estimate for 1985 is 
90 percent. The software costs for the U.S. manned space 
missions from 1960 to 1970 have been estimated at 1 billion 
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dollars. Ref. 1. In the private sector the same relative 
magnitudes and trends are present. The development cost of 
a modern general purpose computer is about one half for 
hardware and one half for software (operating systems, 
compilers, support programs). The operating system 
developed by IBM for the 360 series is estimated to have 
cost 200 million dollars or about 1000 dollars per 
instruction in the final version. 

The production of software can be divided into three 
phases: 

* analysis and design, 

* writing programs and 

* test and integration. 

Data on how time, effort and money are divided among these 
three phases gives some indication of why software 
production is so costly. The fraction of time, effort and 
money for each phase differs from application to 
application; however, data from some large projects show 
similar experience. Estimates are given in Refs. 1 and 2 for 
some military command and control systems: analysis and 
design is about 35 percent, writing programs 15 percent and 
test and integration 50 percent. For space projects the 
estimates are 35, 20, 45 percent. For the IBM 360 operating 
system the estimates are 35, 1 5,. 50 percent. Data for 

business applications indicates less for testing and 
integration and more for analysis and design than the above 
data. The surprising amount of time, effort ani money for 
test and integration is often the item most underestimated 
in planning computer projects, as described in Ref. 3. 
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B. DEFINITIONS 


In the field of software engineering there is little 
agreement on the definition of terms, such as the definition 
of software reliability. In order to make the understanding 
of this paper easier the following definitions will be 
adhered to as much as possible. 

1. Te r ms 

So ftware re li ab il ity is the probability that a 
computer program will perform its intended function for a 
specified interval under the stated conditions according to 
Ref. 4. 


Te stin g is an effort to determine the presence of 
software errors, not their absence. 

Sof tware error is an error in programming logic 
which leads to undesirable results during program execution. 

Module is a particular physical combination of 
program instructions that is independent of others with 
respect to compiling, assembling and loading. 

Program is a set of modules. 

Progra m complexity may be described by 
characteristics such as program size, incidence of branch 
instructions, incidence of loops, incidence of subroutine 
calls and variety of instructions. 
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No n -br anch instructions may be either 

computational or input/output instructions. 

Struct ured Erogramming is a programming technique. 
Ref. 5, in which a program with one entry and one exit can 
be written using only the following programming 
progressions: 

* Sequence 

* IF THEN ELSE 

* DO WHILE 

Direct ed araph is a geometric graph, consisting of 
nodes and arcs, with a direction of traversal associated 
with each arc. 


C. CLASSIFICATION OF ERRORS 


Software errors are classified as follows: 

* Mistakes in logic at the flow chart level, 

* Computation and assignment, 

* Sequencing and control, 

* Input/output, 

* Declarations, 

* Keypunching/clerical errors committed in writing 

instructions on coding sheets, 

* New errors introduced as a result of design changes: 

■ unexpected side effects caused by changes, 

■ logical flaws in change to design, 

■ inconsistencies between changed design and 

implementation, 

■ inconsistencies in original and changed hardware 
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TESTING AND ERROR DETECTION 


The life cycle of a program is composed of the following 
phases: 

* Design and analysis, 

* Module development and testing, 

* System integration testing, 

* Functional testing, 

* Maintenance. 

The cost of error detection and repair during system 
integration testing is three times that of testing an 
individual module during module development testing at TRW, 
according to Ref. 6. Therefore, the objective should be to 
reduce the number of errors detected during system 

integration testing and increase the number (proportion) 
discovered during module development testing. 

In many moderate and large computer projects, a 

programmer writes and debugs a module and then gives it to a 
test group. The test group tests the module, integrates it 
with other modules and then continues testing. The module 
is tested by supplying an input to the module and then 
comparing the outcome to the known correct outcome. If 
there is a mismatch between observed and correct output, an 
error has been detected. When an error is detected the 
module is given to a programmer who locates and corrects the 
error and then returns the module to the test group. Notice 
the distinction between testing, which is supplying inputs 
and observing outputs, and debugging, which is the highly 
individualized detective work needed to locate and correct 
errors. In debugging, the programmer needs a detailed 
knowledge of the structure and operation of the module. The 
tester is frequently unaware of module structure and 
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operation; he needs only to understand the function of the 
module. 

Most computer programs have a large number of potential 
inputs; each may exercise a program in a different way. The 
sequence of instructions of the program that results from a 
particular input is called the "path" or "thread" associated 
with that input. Testing by submitting inputs to the 
program checks only the paths associated with those inputs. 
For programs with a very large number of inputs, testing can 
be only a relatively small sampling of all possible inputs, 
as described in Ref. 3. 
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ERROR DETECTION MODEL 


A. NEED FOR A MODEL 


Testing is a critical part of software projects because 
it measures and affects the final quality of the software 
and it consumes a large part of project time and resources. 
Testing also reveals the strengths and weaknesses of the 
analysis, design and coding of the software and gives an 
estimate of the success or failure of the software after 
release. Thus it is important to understand the testing 
process and to understand the relationships between testing 
and the various decision variables that may be controlled 
during analysis, design and coding. 

A difficult facet of program testing involves the 
selection of inputs. The tester, who generally is not the 
person who wrote the code, does not know the specific path 
that an input will execute. Presently there is no software 
tool, such as an automatic test data generator, that would 
allow the tester to force an input to follow a certain path. 
The closest system designed for testing a certain path is 
interactive, where the tester selects whichever instruction 
is to follow the previous one. In this way a particular 
path is followed. Ref. 7. This is obvioulsy a slow and 
cumbersome way to check out all, or many, of the possible 
paths in a program. 
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Obviously, inputs should be chosen so that a high 
percentage of the critical paths of the program will be 
exposed to testing. However, this objective must be weighed 
against the cost of machine time for debugging and the cost 
of programming personnel for error correction. A related 
matter is the determination of when to stop testing. It is 
usually infeasible to subject a program to all possible 
input combinations because of resource constraints. Various 
software packages are available for recording and analyzing 
the following types of data: count and frequency 
distribution of types of instructions executed; indication 
of code which is not executed; and indication of code which 
is impossible to reach. Ref. 8. Although this type of 
instrumentation is helpful for tracing program behavior, 
once a set of inputs is selected, it does not solve the 
problem of selecting the number and type of inputs in the 
first place. 

Thus, there is a need for a model to examine the 
relationships between the number of inputs and paths 
traversed, for a given program structure, and the number of 
remaining errors, fraction of the program exposed to 
testing, execution time and repair time. It is of interest 
to determine the number of inputs required to achieve a 
specified number of remaining errors for various structures, 
when the same number of original errors is used with each 
structure. In addition it is desirable to identify 
programming structures which have complexities that make it 
difficult to detect errors. 
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B. BASIC MODEL DESCRIPTION 


% 


» 


1. Model Characteristics 


Program complexity may be described by 
characteristics such as program size, incidence of branch 
instructions, incidence of loops, incidence of subroutine 
calls and variety of instructions. Another view of program 
complexity can be obtained by considering the structure of 
the program to be a series of nodes, arcs and loops in the 
form of a directed graph as shown in Figure 1. 


In the directed graph used in the simulati 
nodes represent connection points where parts of th 
may merge and/or branch and arcs represent a se 
nonbranching instructions such as computati 
input/output. Instructions are located in arcs a 
are located in some of the instructions. An input 
path from the start node to an exit node. Beginni 
start node an input causes execution of the instruc 
its path, consuming test time, until, an 
encountered. After the error is thus detected 
repaired, consuming repair time. There is, hove 
risk that the repair will introduce a new error 
instruction. Restarting at the initial node, exe 
begun again with the same input. This process is 
until there are no errors on the path. 


on model, 
e program 
guence of 
on and 
nd errors 
defines a 
ng at the 
tions on 
error is 
, it is 
ver, some 
in some 
cution is 
repeated 


Some relative measures of program complexity which 
are applicable to a directed graph representation of program 
structure are: 
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Fig 1 Directed Graph Representation of 


Program 
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* ratio of actual number of arcs to the maximum 

possible number of arcs (every node is connected 
to every other node by an arc), 

* ratio of nodes to arcs, 

* ratio of loops to total number of arcs. 

The size of a program is a measure of complexity in an 
absolute sense. In terms of a directed graph structure, 
size is determined by the number of nodes, which establishes 
the number of branch points in a program, and by the number 
of arcs, which establishes the degree of straightline coding 
between branch points. Increasing values of the above 
relative and absolute measures represent increasing program 
complexity. 

2. Model Simu l at ion 

The error detection model was written in FORTRAN IV 
and has been developed and used on the Naval Postgraduate 
School’s IBM 360/67 computer. The program has been executed 
40 times in the production mode. The simulation program 
consists of 639 FORTRAN statements, requiring 194,000 bytes 
of main memory and executes in 40 to 55 seconds, depending 
on the type of simulation involved. The directions for use 
of the error simulation program are listed in Appendix B, 
Appendix C is the flow chart, which was produced by the 
FLOWCH routine which is part of the Programming Aids Library 
at the Naval Postraduate School, and the FORTRAN code for 
the simulation model is in Appendix D. The directed graph 
was input to the simulation as a node-arc incidence matrix. 
Lacking detailed information about the distributions of the 
pertinent variables in actual systems, there were no 
statistical dependencies among the variables established. 
Thus the random variables were chosen to be independent and 
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to possess the Markov property. This also makes the model 
more tractable for obtaining an analytical solution. 

The number of instructions per arc is an independent 
exponential random variable truncated to an integer. Errors 
are inserted by making the number of instructions between 
errors an independent exponential random variable, which 
results in a Poisson distribution of errors per interval of 
instructions. Errors are inserted by scanning the arcs of 
the node-arc incidence matrix by columns until the count of 
instructions from the last error equals the random number. 

An input is a sequence of random numbers that 
determine which arc to traverse at each branch node. For 
each branch node the probability of taking each arc is 
equal. This could be changed to test the sensitivity of 
error detection to different branch probabilities. 

The repair times for errors are exponentially 
distributed. If many programmers work on error repair with 
each repairing only a small number of errors, the effect of 
experience on error repair may be small so that a constant 
repair rate corresponding to the exponential distribution 
would be appropriate. If few programmers work on repairs, 
experience would be a factor and an increasing repair rate 
distribution would be appropriate. For example the 
log-normal is sometimes used to represent the distribution 
of hardware repair times. Ref. 8. 

The execution times of instructions are 
exponentially distributed. It was assumed that the 
execution time of an instruction does not depend on past 
instruction times. This assumption may not hold if the 
programmer tends to sequence his instructions in certain 
patterns. 
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When errors are repaired, the potential introduction 
of new errors is simulated. New error insertion is based on 
the ratio of the number of instructions changed by error 
repair to the total number of instructions in the arc. The 
arc where the new error is to be inserted is determined on 
an equal probability basis. 


The simulation is 
or its parameter can be 
sensitivity analysis. The 
a significant effect on the 
structure; however, since 
results on a relative basis 
choice of distributions does 


written so that any distribution 
changed for the purpose of 
choice of distributions may have 
simulation results for a given 
the objective is to evaluate 
across various structures, the 
not seem to be critical. 


For each input, data is collected on the number and 
location of errors detected, number and location of new 
errors, number and location of remaining errors, number of 
arcs traversed, time to execute instructions and time to 
repair errors. 


The simulation model was written so that it 
possible to generate random times for each in 
executed and for each error repaired as the s 
proceeds. However, if the instruction times an 
times are independent and identically distri 
decribed above, then it is possible and comput. 
desirable to count the number of instructions exe 
the number of errors repaired and multiply these 
average instruction executing time and the aver 
repair time, respectively, in order to obtain a v 
estimate of each total time. 
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C. MODEL ASSUMPTIONS 

A basic assumption of the model is that the tester has 
some knowledge of the program structure, but that for a 
given input he does not know the specific path that it will 
execute. In actual software projects the test group has 
flow charts and program listings; however, it is infeasible 
to analyze this information because it may contain, thousands 
of lines of coding. Because of the size of the program, the 
complicated internal logic and the large number of paths, 
the relationship between inputs and outcomes is rarely 
understood. One example is in the testing and maintenance 
of large operating systems. The relationship of inputs to 
outcomes is so poorly understood~that even after an error 
has been detected it is often difficult to determine an 
input that will reproduce the error. 

A further assumption of the model is that the tester 
gains no information as the testing proceeds that will 
influence his choice of subsequent inputs. In actual 
software projects the tester should try to make best use of 
any information gained during testing. Various software 
packages are available for recording the following types of 
data count and frequency distributions of instructions 
executed, indication of code that is not executed and 
indication of code that is impossible to reach, Ref. 9. 
However, there are other factors that may make it difficult 
to effectively use the information gained during testing. 
For example, the test plan may be specified in advance with 
no modifications allowed or inputs may be restricted to 
those that will be typical for the program in actual 
operation. For these reasons the model assumptions seem 
reasonable as applied to functional testing. 
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The probability distributions which were used are listed 
below. 

P roperty or Event_Probability Distribu ti on 


* Instructions per arc 

* Instruction execution time 

* Original error occurrence 

* Time to repair an error 


Exponential 

Exponential 

Exponential 

Exponential 


Number of instructions 

Uniform 

affected by repair 

New error occurrence 

Uniform 

(based on 

instructio 

changed/ 

instructio 

Iterations per loop 

Uniform 

Arc selected for new 

Uniform 

error insertion 

Arc selected at branch 

Uniform 


point for traversal 


Since little is known about the type of probability 
distribution which is associated with the above program 
properties and execution events, the selection of 
distributions was, of necessity, based on assumptions. 
However, it was felt that the assumptions were reasonable. 
For example, the seeding of original errors was based on the 
number of instructions between errors being exponentially 
distributed, or equivalently, the presence of an error was 
independent of the presence of other errors. A second 
example was that instructions were placed in arcs according 
to an exponential distribution, or equivalently, the number 
of instructions between branch points was exponentially 
distributed. This implies that the number of instructions 
between two branch points was independent of the number of 
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instructions between other branch points. Although the 
choice of distribution may have a significant effect on the 
simulation results for a given structure, the objective was 
to evaluate results on a relative basis across the various 
structures so that choice of distribution was not critical. 
Although it was possible to vary both the type of 
distribution and its parameters, the usual procedure was to 
keep these factors constant and vary program structure, 
number of inputs and input traversals. 

D. MODEL USES 

The model can be used to influence software design 
decisions by making it possible to compare the error 
detection characteristics of alternative program structures. 
This is valuable, since error detection characteristics are 
good indicators of the time and resources consumed by 
testing. The design flow charts and estimates of branch 
probabilities and number of instructions can be used to 
specify programs in the form of a directed graph. The 
program is then seeded with errors and subjected to random 
inputs. 

The model can also be used to identify the measure or 
measures of complexity that best predict the ability to 
detect errors. To do this it is necessary to gather data 
from the model on the error detection characteristics of a 
variety of different structures and then do a statistical 
analysis. This would make it possible to measure the 
complexity of different programs and then compare the 
estimates of error detection characteristics. Although some 
data has been generated, further work is necessary to 
identify good measures of complexity. 

/ 
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There are other situati 
to compare structures, 
the cost of adding some add 
The results of the mode 
detection characteristics 
structure. The problem of 
structures of different siz 
addressed. 


ons where it is useful to be able 
A frequent problem is to evaluate 
itional feature to the program. 
1 can be used to compare error 
of the original and modified 
how to allocate test effort among 
e and complexity can also be 
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IV 


ANALYSIS OF SIMULATION RESULTS 


A. THE EFFECT OF INCREASING THE NUMBER OF INPUTS 


1. Mod el T est ing 


One would expect that initially there are many 
errors detected in a program with each input and then the 
number of errors detected decreases as additional test 
inputs are used, because much of the program is exposed to 
testing initially. This is illustrated in Figure 2. The 
percent residual errors decreased stepwise as the number of 
inputs increased. In the testing of actual software, after 
finding many errors, there may be long periods of time with 
no error detection followed by a new group of detected 
errors. 


Recall that each input in the model detected all the 
errors in its path, from the input node to one of the output 
(terminal) nodes. In order to explain the stepwise action 
in Figure 2, it must be realized that although the paths 
through the program were, in general, different from 
previous paths, portions of these paths may have involved 
only arcs that had been traversed previously. The model had 
well defined steps where no new arcs were tested for a 
number of unique input paths, as shown in Figure 3. Thus 






30 nodes, 50 arcs, 6 loops 
18 original errors, 11 added errors 



Fig 2 The Effect of Increasing Inputs on Residual Errors 

it can be seen how a new group of errors was detected when 
the model tested previously untested parts of the program. 
Just because an arc has been previously tested does not 
imply that it was error free. As each new detected error 
was repaired, there was soma small probability that a new 
error was introduced in some other portion of the program. 
This newly inserted error may have been inserted in a 
previously tested arc. A check was made on the coverage of 
the arcs by the simulation model. The 
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30 nodes, 45 arcs, 6 loops 
28 original errors, 17 added errors 
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Fig 3 The Effect of Increasing Inputs on fires Tested 

structure checked had 30 nodes, 40 arcs and 6 loops as shown 
in Figure 4. The numbers along the arcs indicate the number 
of times the arc was traversed. For example, the source arc 
at the top was traversed 50 times, or there were 50 
different inputs. Every time an input reached a node it had 
an equally likely opportunity to select any one of the arcs 
emanating from the node. The simulation results bear this 
out as Figure 4 illustrates, where the 50 inputs traversed 
the top arc and split below with 25 going to the left and 25 
going to the right. The arcs which had a backward pointing 
arc or loop around them were traversed more times as shown 
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Fig 5 FORTRAN Program Directed Graph 
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by the number alongside the arc representing the sum of the 
number of times the backward loop was selected and the 
number of times the input arc was selected. Looking at the 
far right hand loop, seven inputs came into the node from 
above, twelve inputs came into the node from the loop and 19 
of the inputs exited the node. The parenthesized numbers 
indicate the number and location of errors. The number of 
errors includes the errors initially seeded and the errors 
inserted when repairing detected errors. 


2. S im u lat i on E xam ple on a Real P ro gram 


In order for the simulation model to be of any 
practical use, it had to be tested on a real program. 
Appendix A contains the code and structure of a textbook 
FORTRAN program for computing Bessel Functions. The column 
labeled "node" corresponds to the nodes in the directed 
graph representation of the program in Figure 5. This 
particular program was selected as an example of a good 
computational program, since it was presented in a numerical 
analysis text. Ref. 10, and an example of a poorly coded 
program, since a casual reading of the code showed a lack of 
use of structured programming technigues. Another reason 
for the selection was that the program could be broken down 
into 30 nodes, which was the same as the test structures. 
It also fit within the range of structures tested having 43 
arcs and 9 loops. The first number inside of the 
parenthesis represents the number of instructions in the arc 
and the second number represents the number of errors, 
original errors plus added errors, in the arc. 

Fifty randomly selected inputs were run through the 
structure. With 16 errors initially seeded, five errors, or 
16.7 percent of the total errors seeded, still remained 
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after fifty inputs. Comparing this result with a test 
structure with 45 arcs and 6 loops and another test 
structure with 44 arcs and 10 loops, the percent residual 
errors in the FORTRAN program was high, illustrated by 
Figure 6. By analyzing the paths each input traversed it 
was noted that six of the nine loops in the FORTRAN program, 
all emanating from the bottom of the graph, going to the top 
of the graph, were very seldom used. Thus some individual 
inputs were not given an opportunity to loop back up to the 
top of the graph, and thus test more branches for a given 
input. 

FORTRAN Program 
30 nodes, 43 arcs, 9 loops 
16 original errors, 14 added errors 



Fig 6 Residual Error Pattern 
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Percent Arcs Tested 


This was borne oat by the results in Figure 7, which showed 
that the percentage of arcs tested was lower for the FORTRAN 
program than for the other structures tested. 


FORTRAN Program 
30 nodes, 43 arcs, 9 loops 
16 original errors, 14 added errors 



FORTRAN 


30 

Number of Inputs 




50 


Fig 7 Arcs Tested Pattern 

B. THE EFFECT OF INCREASING THE NUMBER OF ARCS 

Intuitively, given two programs with the same number of 
nodes, and a different 'number of arcs emanating from the 
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nodes, one would expect that the program with the greater 
number of arcs, or the more complex program, would have the 
higher percentage of residual errors. By the same reasoning 
one would expect the more complex program to have fewer arcs 
tested with a given number of inputs. 

Fifty random inputs were used on each of the program 
structures described below. Each structure contained thirty 
nodes and six loops. Retaining the concept that each node 
represents a branch or decision point in the program, the 
most simple structure that can be defined, using thirty 
nodes, must have a minimum of forty arcs. By definition, to 
establish a node there must be at least three arcs, in any 
combination, either terminating or emanating from the node. 
Thus, the minimum number of arcs in a structure is 3/2(the 
total number of nodes minus the number of entry and terminal 
nodes). Recall that an arc was defined as either a forward 
or backward pointing arc, called a loop. Starting with 
forty arcs and adding five more to each structure, five 
structures were simulated with 40, 45, 50, 55 and 60 arcs. 
After fifty inputs the percent residual errors increased as 
the number of arcs increased, as Figure 8 illustrates. The 
percent residual errors was chosen as the vertical axis in 
Figure 8 rather than residual errors since the number of 
errors, original errors plus added errors, varied in each of 
the five structures. The reason for the variable number of 
errors was that each time a new structure was defined, the 
error simulation program would randomly seed all the 
original errors again, thus errors could have been inserted 
into the added arcs. 
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30 nodes, 40 


60 arcs, 6 loops 



Fig 8 Relation Between Complexity and Residual Errors 

Similarly, Figure 9 illustrates the effect of increased 
complexity on the percentage of the arcs tested. As the 
number of arcs increased, the percent of the arcs tested 
decreased. 

Examining the paths traversed by each input gave some 
insight as to why an increased number of arcs caused higher 
residual error and lower percentages for arcs tested. When 
an arc was added, the number of arcs emanating from a node 
increased. There was a positive probability that the added 
arc could contain an error as the entire structure was 
seeded with errors anew. As the number of arcs increased, 
there were also more arcs which provided shorter paths to an 
exit node 
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30 nodes, 40 - 60 arcs, 6 loops 



Fig 9 Relation Between Complexity and Percent Ares Tested 

by connecting a node closer to the input with a node closer 
to one of the outputs, thus leaving some intermediate arcs 
untested. 

Repair time turned out to be unrelated to complexity. 
The number of errors initially seeded controlled the repair 
time. These results can be seen in Figure 10, where the 
data points have been smoothed to give straight lines. 
Several structures are shown in the plot of repair time 
versus percent residual errors. The amount of time required 
to repair errors, for a given percentage of residual errors, 
increased as the number of errors initially seeded 
increased. 
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40 - 60 arcs, 6 loops 


30 nodes. 



Fig 10 The Effect of Arcs on Repair Time 


Generally 
tested and the 
approximately 
increased, the 
decreased. In 
band of values, 
given percenta 
mean value. 


the relationship between the percent arcs 
percent residual errors can be described as 
linear. As the percentage of the arcs tested 
percentage of the residual errors remaining 
Figure 11, the shaded area represents the 
corresponding to various structures, for a 
ge of arcs tested. The solid curve is the 
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30 nodes, 40 - 60 arcs, 6 loops 



Fig 11 The effect of Complexity on the Relationship between 
the Residual Errors and Arcs tested 

C. THE EFFECT OF INCREASING THE NUMBER OF LOOPS 

Improper loop indexing is usually near the top of a list 
of most frequently occurring errors. Ref. 1. Many people 
think that loops should be eliminated as a program 
structure. Note that the only influence that loops play in 
this model is with respect to coverage. The model does not 
account for errors in the loop counter or failure to get out 
of a loop. One of the results of the analysis, as shown in 
Figure 12, was that an increase in the number of loops had 
no significant effect on the percentage of residual errors. 
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14 and 20 


Structures with zero, 5, 6, 10, 14 and 20 1 

analyzed using the error simulation program. r 
residual errors was chosen as the vertical axis ra 
residual errors since the number of errors, origi 
plus added errors, varied in each of the six s 
The reason for the variable number of errors was 
time a new structure was defined the error 
program would randomly seed all the original err 
thus errors could have been inserted into the adde 
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30 nodes, 34 


54 arcs, 0-20 loops 



Fig 12 The Effect of Loops on Residual Errors 

The reason for the independence of the percent residual 
errors from loops can probably be explained by the fact that 
there was no distinguishaole difference between the percent; 
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of arcs tested in the six cases with 0, 5, 6, 10, 14 and 20 
loops. This is shown in Figure 13. 

30 nodes, 34 - 54 arcs, 0-20 loops 



Fig 13 The Effect of Loops on Percent of Arcs Tested 

By examining the paths the inputs trace, the explanation 
of the above becomes obvious. After an input completes a 
loop, it once again has an opportunity of branching out of 
the loop, thus testing more arcs than a structure with no 
loops. Each time another loop was added, the probability of 
branching out of all t*he loops increased at approximately 
the same rate as the increased number of loops. This 
concept was reinforced by the data shown in Figure 14. 
After the structure was expanded to nine loops, additional 
loops had no effect on error detection. The percent of the 
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total residual errors in the structure that resided in the 
loops was a constant 59 percent and the percent in the arcs 
was a constant 41 percent for structures with nine to twenty 
loops. This data was derived by starting with a structure 
containing 20 loops, seeding errors, and then analyzing the 
structure. The error simulation program would then delete a 
loop and its associated error, if one had been seeded, and 
then analyze the new structure. 


Error locations constant 
Variable number of errors 



Fig 14 Breakdown of Residual Errors in Loops and in Arcs 

It required essentially the same amount of repair time 
to decrease the percent residual errors to a certain level 
for all the structures containing loops, as shown in 
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Figure 15. The structure with no loops took less than half 
the repair time to get to the same percent of residual 
errors as a structure with loops, for a percent residual 
errors less than 50. 


30 nodes, 34 - 54 arcs, 0-20 loops 
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Fig 15 The Effect of Loops on Repair Time 

In Figure 10 it was shown that the number of errors 
initially seeded determined repair time. In Figure 15, the 
structure with no loops had fewer initially seeded errors 
than all the structures with loops, so the above 
relationship holds for this case. However, the number of 
errors initially seeded had no distinguishable effect on the 
repair time of the structures with loops. Thus, the repair 
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time must be a function of both the number of errors 
initially seeded and whether or not loops are present. 

It was not possible to make any judgements concerning 
the determinants of execution time. This was due to the 
fact that all but one structure tested had loops. Loops 
were executed a variable number of times as determined by a 
uniform distribution which established the number of 
iterations. The effect of a doubly nested DO loop was 
captured by allowing an input to have an equiprobable chance 
of branching back up to the start of the loop or of 
branching farther down the structure. The relationship 


30 nodes, 34 arcs, no loops 
18 original errors, 6 added errors 
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Fig 16 Execution Tine for a Structure with No Loops 
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between the number of inputs and the cumulative execution 
time for a structure with no loops is examined in Figure 16. 
A plot for any of the structures with loops has points 
scattered all over due to the random effect of the loops on 
execution time. 

The effect of the percentage arcs tested, with loops 
present, on the percent residual errors can loosely be 
described as a linear relationship. As long as the 
structures all had loops, the curves of the percentage of 


30 nodes, 34 - 54 arcs, 0-20 loops 



Fig 17 The Effect of Loops on the Relationship between 
Residual Errors and Arcs Tested 
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arcs tested versus the percent residual errors all fall 


within a narrow band of values as shown by the shaded area 
of Figure 17. The curve for a structure with no loops is 
also plotted. 

The mean of the percent residual errors of structures 
with a variable number of arcs and a constant number of 
loops and the mean of percent residual errors of structures 
with a variable number of loops are plotted in Figure 13, 
which shows that the two curves are almost identical. 
However, as Figure 18 illustrates, the structure with no 
loops required significantly fewer arcs to be tested to 
achieve the same level of residual error percentage as 
compared to the structures with loops. 
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Fig 18 The Advantage of No Loops 
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D. 


REPLICATING A SINGLE INPUT 


1. Model T es ting 

The usefulness of the model is now examined for 
predicting the ability of detecting errors in an actual 
program. Four pieces of information are of importance for a 
manager conducting module development testing of computer 
programs. These are: the percent or number of residual 
errors, the percent or number of arcs tested, the amount of 
repair time, and the amount of execution time. 

The simulation model was used on ten different test 
structures to see if this information could be predicted. 
For each structure a single randomly selected input was run 
and the above data was collected. This process was 
replicated 100 times, or in other words, 100 randomly 
selected inputs were used with each input using the same 
structure and the same number of errors seeded in the same 
places. Statistics such as mean, median, variance, standard 
deviation, etc., were calculated. 

As an example, the basic 30 node, 40 arc, 6 loop 
structure had 24 errors initially seeded. The simulation 
model produced a mean of 78.79 percent residual errors with 
a standard deviation of 9.10 for one input. Thus, one could 
estimate that based on 24 original errors, 78.79 percent of 
the errors will remain after one input. Similar statistics 
were determined for percentage of arcs tested and repair 
time. Execution time was found to have a high variance. 
For instance, the mean execution time for one input for the 
above structure was 32.50 seconds with a standard deviation 
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of 50.15. Thus, estimates of execution time based on the 
mean would be subject to high error. 


2• Si m ulation Exam ple on a Real Program 


The simulation moiel was used on the FORTRAN Bessel 
Function program described earlier. It was found that, 
based on 16 original errors the expected percent residual 
errors was 84.26 percent with a standard deviation of 9.09, 
or 15.74 percent of the original errors could be expected to 
be found and corrected with one input. Similarly, 17.70 
percent, with a standard deviation of 8.65, of the arcs 
could be expected to be tested by one input. Of prime 
importance to the project manager, 1.41 hours of repair 
time, with a standard deviation of 1.18, a relative measure 
for the manager to use when comparing alternative 

structures, could be expected to be devoted to detecting and 
repairing 15.74 percent of the errors. 

E. THE EFFECT OF COMPLEXITY 

The following complexity measures will be used: 

* AMA is the ratio of the number of arcs in the 

structure to the maximum number of arcs possible for the 
given number of nodes, 

* NA is the ratio of the number of nodes to the number 
of arcs in the structure, 

* LA is the ratio of the number of loops to the number 
of arcs in the structure. 

Using these relative complexity measures, it was of 

interest to see how each of the measures affected the 

percent residual errors and the percentage of arcs tested. 
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Five different structures with a constant number of loops 
and a varying number of arcs and six different structures 
with a varying number of loops were examined. For each 
structure, 100 replications of a single input were simulated 
using the error simulation program, and statistics were 
gathered about the percent residual errors and the 
percentage of arcs tested. 

In Figure 19, the percent residual errors after one 
input increased as the complexity increased. In this case 
the complexity measure was the ratio of the actual number of 
arcs to the maximum number of arcs possible with a given 
number of nodes. Similarly, using the same complexity 
measure the percent arcs tested after one input decreased as 
the complexity increased, as shown in Figure 20. In both 
Figures 19 and 20, the standard deviation from the mean. 



Fig 19 The Effect of AMA on Residual Errors 
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0 


represented by the dashed lines# decreased as complexity 
increased. 



Ratio of Nodes to Arcs 


Fig 20 The Effect of AI1A on Arcs Tested 

Using the ratio of nodes to arcs as a complexity measure 
similar results were obtained. The percent residual errors 
increased and the percentage of the arcs tested decreased as 
the complexity increased. These results can be seen in 
Figures 21 and 22# where increasing complexity is from right 
to left. Note that there was an even sharper decrease in 
the standard deviation as complexity increased in both 
Figures 21 and 22. 
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Variable number of arcs 



Fig 21 


The Effect of N A on Residual Errors 
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Constant number of nodes 
Variable number of arcs 


+> 

I 

M 


n 

V 

-p 

•O 

a> 

-p 

to 

a) 

H 

<0 




Fig 22 The Effect of NA on Arcs Tested 

A third complexity measure used was the ratio of loops 
to arcs. In Figure 23, the percent residual errors 
remained constant for one input as the complexity increased. 
Once again this reinforced the idea that loops expose more 
arcs to testing at the same rate as the additional arcs 
increase the complexity. In Figure 24, it can be seen that 
the added complexity had no effect on the percentage of the 
arcs tested after one input, with the mean after 100 
replications being a constant 22 percent. Since loops were 
also defined as arcs, the ratio of loops to arcs did not 
increase linearly as the number of loops increased. 


* 
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Percent Arcs Tested after 1 Input 



Fig 23 The Effect of LA on Residual Errors 
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Percent Arcs Tested after 1 Input 



Fig 24 The Effect of LA on Arcs Tested 












V. 


CONCLUSIONS 


A. INFERENCES THAT MAY BE DRAWN FROM THE RESULTS 


This error simulation model has the potential to be 
developed into an aid to a manager for the allocation of 
test effort among the modules of a program. Since it is 
impossible to test all possible paths in even a moderate 
sized program, due to the large number of possible paths 
created as the number of node's (branch points) increases, it 
is desirable to know which areas of the program should be 
tested more thoroughly than others. By writing the program 
as a directed graph, a manager may use the simulation model 
to give an estimate of how long to test each module, by 

using the relationship between error detection and number of 
inputs. 

The simulation model has shown that a structure with 
many branches is more likely to have a higher percentage of 
residual errors than a structure with straight line code. 
In certain instances, a structure with many loops is 
preferable to one with many branches. The model verified 
this conclusion in terms of percentage of the arcs tested, 
where the added loops were an advantage in testing more 
code, whereas additional branches left many arcs untested. 
Thus, this model will identify structure types for the 
manager to select for more thorough testing. 
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Another feature of the model is the ability to estimate 
residual errors, percentage of arcs tested, and the amount 
of repair time required by replicating the testing of a 
structure with different inputs. By using this type of 
simulation for the structures which the manager must test, 
he would obtain useful information for allocating the test 
effort. 


B. FUTURE WORK 


The 
phase 
module 
prior 
of the 
used f 
defined 
and the 


error simulation model should be most useful in the 
of testing computer programs known as individual 
development testing. Each module should be tested 
to systems integration testing. With a modification 
definitions presented herein, the model could be also 
or systems integration testing, if the modules are 
as the nodes of the directed graph and the inputs to 
outputs from each module are defined as the arcs. 


One use of the simulation model would be to write a 
program in the "normal" fashion and then write a program 
using structured programming techniques. The two structures 
could be simulated with the same numbers of initially seeded 
errors and inputs, and then comparisons of the percent 
residual errors, percent arcs tested and repair time could 
be made. 
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APPENDIX A 


FORTRAN EXAMPLE PROGRAM CODE 


N od e 


1 


2 


4 

7,11 

8 

12 

15 

14 

17 

21 


C 

C 

c 


c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Code 

"3TJEIT0UTINE FOR COMPUTING J-N, Y-N, I-II- AND K-N 
ASYMPTOTIC EXPANSION IS USED FOR LARGE ARGUMENT 

SUBROUTINE BESSEL (X, NORD,BJN,BY,BIN,BK) 

PI = 3. 14152927 ' 

GAM = .57721566 

FN = NORD 

IF (X - FN - 6.) 1,200,200 

1 XA = X / 2. 

XB = XA * XA 

COMPUTATION OF J-ZERO AND I-ZERO 
N = 0 


COMPUTATION OF J-N AND I-N BY POWER SERIES 


3 

10 

12 

20 


25 

30 


9999,20,12 
XA / AN 

- 1 . 


AN = N 
T = 1. 

S = -1 , 

IF (AN) 

T = T * 

AN = AN 
GO TO 10 
BJN = T 
BIN = T 

DO 30 K = 1,K 
DEN = K * (K + N) 

T = T * X3 / DEN 
IF ( (3JN + T) - BJN) 
BJN = BJN + T*3 
BIN = BIN + T 
S = -S 


25,50,25 


50 IF (N - 1) 75,130,55 

K t- N IS COMPUTED FROM ASYMPTOTIC EXPANSION 
IF X IS GREATER THAN N + 3 

55 IF (X - FN - 3.) 1111,1111,200 

CALCULATION OF J-1 AND 1-1 

65 N = 1 

BJO = BJN 
BIO = BIN 
BYO = BY 
BKO = BK 
GO TO 3 
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Node 

r* 


Code 


c 


COMPUTATION OF K-ZERO AND Y-ZERO 


c 

75 

BY - 2. /PI * (GAM + LOGF(XA)} * BJN 

BK = - (GAM + LOGF (X A) ) * Bx N 




T = XB 

S = 1. 

XI = 1. 

22 



DO 110 K = 2,K 

AK = K 

25 


1 00 

IF ((BY + T*XI) - BY) 100, 120, 100 

BK = BK + T*XI 




BY = BY + 2./PI*S*T*XI 

T = T*XB / (AK*AK) 

XI = XI + 1./AK 

28 

c 

110 

s = -s 

27 

c 




c 



29 


120 

IF (NORD) 9999,55,65 


c 


COMPUTATION OF Y-1 AND K-1 BY WRONSKIAN 


c 

130 

BY = (BJN*BYO - 2./(PI*X)) / BJO 

BK = (1. / X - BIN*BKO) / BIO 


c 



c 


Y-N AND K-N BY RECURSION FORMULAS 


c 


FOR ORDERS HIGHER THAN ONE 


c 


P = 1. 

20 


140 

IF (NORD - 1) 9999,55,140 

BY 1 = BY 

BK1 = BK 

BY = 2„*P / X*BY1 - BYO 

BK = 2.*P / X*BY1 + BYO 


c 


P = P + 1. 


c 



23 


146 

I? (NORD - 2) 9999,150,146 

BYO = BY1 

BK = BK1 

NORD = NORD - 1 

GO TO 140 


c 

150 

N = P 

GO TO 3 


FORMULAS 
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Code 


Node 


3 

5 

10 

12 

13 


16 


18,24 


6,9,19, 

26,30 


C 

C 

c 

c 

200 


c 


210 


220 

230 

240 

241 
C 

C 

c 


250 


C 

1111 

C 

9999 


COMPUTATION OF J-N, I-N, K-N, AND Y~N 
BY ASYMPTOTIC EXPANSIONS 

C = 4 * NORD * NORD 
D = 8. * X 

CONI = SQRTF (2./(PI*X) ) 

CON2 = 1./SQPTF(2.*PT*X) 

CON3 = SQRTF(PI/(2.*X)) 

AN = NORD 

PHI = X - (2.*AN + 1.) / 4. * PI 
M = X + 1. + SQRTF (X * X + AN * AN) 

T = (C - 1.) / D 

S = 1. 

0 = 1 . 

PN = 1 . 

QN = T 
BK = 1. + T 
BI = 1. - T 

DO 240 I = 2,M 
AI = 1 

T = (C - (2. * AI - 1.)**2) / D*T / AI 

BK = BK + T 

BI = BI + T*S 

IF (S) 220,9999,210 

PN = PN - T*S*U 

U = -U 

GO TO 230 

IF ((QN + T) - QN) 230,241,230 
QN = ON - T*S*U 
S = -S 
CONTINUE 

BK = EXPF(-X) *CON3*BK 

ASYMPTOTIC EXPANSION IS USED ONLY FOR K-N 
IF X IS BETWEEN N + 3 AND N + 6 

IF (X - FN -6.) 1111,250,250 

BJN = CONI * fPN*COSF (PHI) - QN*SINF (PHI)) 
BY = CONI * (PN*COSF (PHI) + QN*SINF (PHI) ) 
BIN = EXPF(X) * CON2 * 31 

NORD = FN 
RETURN 

STOP 


END 
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APPENDIX B 


DIRECTIONS FOR USE OF THE ERROR SIMULATION PROGRAM 


Random numbers are drawn using the Naval Postgraduate 
School Random Number Generator Package LLRANDOM 
(NPS55LW73061 A) . These calls are RANDOM for a uniform 
distribution and EXPON for an exponential distribution. 

Printer plots are drawn by the PLOTP routine which is 
part of the IBM supplied routines. Directed graph plots are 
drawn by the DRAWP routine available on the Naval 
Postgraduate School IBM 360/67 System connected to a CALCOMP 
plotter. The subroutine GRAPHO is limited because it 
requires the directed graph to be constructed in such a way 
that the rightmost leg of the tree is the longest leg, which 
means that the rightmost terminal node is the lowest node. 

Histograms are plotted and analyzed by the HISTG routine 
available on the Naval Postgraduate School IBM 360/67. The 
functions SIN, COS and ATAN2 are standard IBM 
FORTRAN-supplied subroutines. 

The analysis of loops in the simulation assumed that 
there were instructions in the arc from i to j and also in 
the backward arc from j to i. This may not be realistic, 
but does prove useful when simulating the actual path. When 
an input reaches a branch point that has a backward arc 
emanating from it, the model simulates a random number of 
iterations using the number of instructions in the backward 
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arc. After completing the simulated looping, it was 
possible to trace the same path or a different path if there 
were intermediate branch points. If the same path was 
selected, there was another chance for the backward arc to 
be selected. This is a common programming structure, where 
a loop is executed for so many iterations, a few parameters 
are changed and the loop is executed again for another or 
the same number of iterations. 


The comment section of the simulation program, titled 
"Directions for Use", describes the use of each data input 
to the program as well as the common values utilized during 
this research. The numbers along the top of the card are 
the column numbers, the middle row of numbers are sample 
data and the bottom row of numbers are either the 
corresponding node for each of the pieces of data above or. 
data item names. The formats used and how each data card 
was designed are shown below. 



col It 
data 


Fig 25 Each item in paragraph A of the directions for use 
is on a separate card with 13 format except item A (12) which 

is A5 format (***♦*) . 
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Fig 26 Adjacency matrix card - 1615 format. 



data 


node # 


Fig 27 Second card of adjacency matrix - based on a graph 

with 30 nodes. 
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Fig 29 Second card of the matrix of arc lenths - based on a 

graph with 30 nodes. 



col # 

data 

id 


Fig 30 The two nodes which determine the arc to be added or 
deleted from the structure - 215 format. 


10 



col # 

data 

id 


Fig 31 The seed for the random number generator - format 

1110. 
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FLOW CHART 


• CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PAGE 1 

I 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccccc 


100 

108 

140 


20 3 

204 

205 


SIMULATION OF ERRORS IN SOFTWARE 


C 

c 

c 


cccccccccccccccc ccc ccccccccccccccccccccccccccccccccccc ccccccccccc 

COMMON/ALLL/N, N0DES(30,30) 

COMMON/ARRAY/X( 30 ,30) , NUMPTS(30), AST OUT (30 ) , AVC.HNG 

COMMON/E RRQ R/ISEED(30,30) , ITER(30,30), MEANER, INST, MEANIT, IZ 

COMMON/OUT/NSEED, INPUT, SVSEED(30,30) 

C OMMON/G EN/ME ANLN, IX 
COMMON/MEW/IW 
COMMON/GRAF/INREAD 

DIMENSION TTIME(50), TREPI50), TLINKI50), TKOUNT(50), TINPUT(50) 
DIMENSION PSEED(50), F(50), FNA(50), FLA(50) 

DIMENSION TSEED(50), TINSTI50), PTEST(50) 

I NT EGE R*4 CHANGE, AVCHNG, TESTED(30,30), TEST, SIMNUM, SVSEED 
INTEGERS DELADD 
REAL*4 MTTR 
FORMAT (13) 

FORMAT (110) 

FORMAT (1615) 


151 

FORMAT ( ' 1’ 

, 13X, 

•THIS 

RUN 

WILL 

VARY 

THE 

NUMBER 


ETWEEN EACH 

NODE. 

' ) 






152 

FORMAT Cl' 

, 13X, 

'THIS 

RUN 

WILL 

VARY 

THE 

NUMBER 


SED.•) 








153 

FORMAT |'l' 

, 13X, 

•THIS 

RUN 

WILL 

VARY 

THE 

NUMBER 


HE NODES.') 








154 

FORMAT Cl' 

, 13X, 

•THIS 

RUN 

WILL 

VARY 

THE 

NUMBER 

155 

FORMAT ('I' 

, 13X, 

• THIS 

RUN 

WILL 

VARY 

THE 

RANDOM 


DETERMINES 

THE IN 

PUT PATH.•) 





156 

FORMAT (*0« 

, 13X, 

•THE 

NUMBER OF 

INSTR 

UCT 

IONS IN 


READ IN AND 

IS MO 

T RANDOM.') 





180 

FORMAT ('O' 

, 13X, 

•ADJACENCY 

MATRIX'//7X, 

30( 14)) 


FORMAT ('0*, 13X, 'NOTE: NONEXISTENCE OF AN ARC IS INDICATED BY A 

ZERO'//) 


FORMAT ('O', 13X, 'THERE WILL BE 1 , 15, ' 

.'//) 

FORMAT (7(15, 5X)) 

(CONTINUED ON PAGE 2) 
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2 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

I 

FORMAT (A5) 


PAGE 


2C6 
20 7 


208 

209 

227 
230 
24 0 

24 5 
250 

25 5 
260 
265 

270 

275 

280 
281 
28 2 

283 

284 

285 

286 

287 

288 
28 9 

290 

291 

292 


('O', 12X, 'INPUT ASSUMPTIONS 


FORMAT ('O', 12X, * 

N TIME PER INSTRUCTION HAS 


AN 


( VARIABLES) :'/./13X, 'EXECUTIO 
EXPONENTIAL DISTRIBUTION WITH A MEAN 
a, 'MEAN TIME TO REPAIR AN ERROR HAS 
AN EXPONENTIAL DISTRIBUTION WITH A MEAN OF', 13,' MINUTES.'//13X, 
•THE NUMBER OF INPUT PATHS BEING INVESTIGATED IS', 13, '.'//13X, 

•ERRORS ARE SEEDED BY AN EXPONENTIAL DISTRIBUTION WITH A MEAN OF', 


OF', 13,' MILLISECONDS.*//13X 


13 


INSTRUCTIONS PER ERROR.•//13X, 


CTIONS CHANGED WHEN AN ERROR IS FOUND IS', 13,'.') 


'THE AVERAGE NUMBER OF 

• - T -X . « _ I 


INSTRU 


FORMAT (•O' , 
DISTRIBUTED 


13X, 'THE NUMBER OF ITERATIONS PER LOOP IS UNIFORMLY 


FROM 1 


TO 1 , 


13, 


) 


FORMAT ('O', 13X, 'THE MEAN NUMBER OF INSTRUCTIONS 


('O', 13X, ' 

DE IS EXPONENTIALLY DISTRIBUTED WITH A MEAN OF', 

FORMAT ('O', 13X, 'SAMPLE INPUT PATH FR0M»//13X, 'NODE TO NODE'//) 
F ORMAT('O', 11X, 14, 4X, 14) 

FORMAT!*0',13X,»FOR INPUT NUMBER',13,' THE EXECUTION TIME IS'- 


BETWEEN EACH NO 

13, ' .* ) 


F 8.2, 

// 13X, 


• SECONDS WITH A TOTAL REPAIR TIME OF', F8.2, • HOURS' 
•AND AN AVERAGE LINK TRAVERSAL TIME OF', F8.2,' SECONDS') 


FORMAT ('O', 13X, 'ERROR FOUND'//) 

FORMAT( r 1 * » 13X, 14, ' SEEDED ERRORS REMAINING'//13X, 'MATRIX OF 
SEEDED ERRORS REMAINING*//7X, 30(14)) 


FORMAT(* ', 

FORMAT('O', 


7 X ■ 


30(A4)) 

15, • *', 


30(14)) 


FORMAT(' 1 *, 11X, 13, 


INPUTS'//13X, 14, • INSTRUCTIGNS'//13X, 


13, ' ERRORS SEEDED*//13X, 13, ' RESIDUAL ERRORS•//13X, F6.2, 


' PERCENT RESIDUAL ERRORS•//13X, F6.2, 
//) 


• PERCENT OF ARCS TESTED* 


FORMAT('O', 9X, F10.5, • SECONDS EXECUTION TIME'//12X, 

F10.6, • HOURS TO REPAIR ERRORS' //12X, F10.6, ' SECONDS AVE 

RAGE ARC TRAVERSAL TIME'//) 

FORMAT('O', 12X, F5.2, ' IS THE RATIO OF ACTUAL TO MAXIMUM NUMBER 
OF ARCS'// 13X, F5.2, ' IS THE RATIO OF NODES TO ARCS»//13X, 

F5.2, • IS THE RATIO OF LOOPS TO ARCS'//) 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


(*0»//13X, 

( ' 0' //I3X , 

( '0'//13X, 

( *0'//13X»- 


•HISTOGRAM FOR EXECUTION TIME') 
•HISTOGRAM FOR REPAIR TIME') 

•HISTOGRAM FOR PERCENT RESIDUAL ERRORS') 
•HISTOGRAM FOR PERCENT ARCS TESTED') 


Cl 


IX) 


FORMAT (•0*//13X,'EXECUTION TIME VS RESIDUAL ERRORS') 

FORMAT ('O'//13X,'EXECUTION TIME VS PERCENT RESIDUAL ERRORS') 

('0 1 //13X, 'REPAIR TIME VS RESIDUAL ERRORS') 

•NUMBER OF INPUTS VS RESIDUAL ERRORS') 

•NUMBER OF INSTRUCTIONS VS RESIDUAL ERRORS') 
•RATIO OF ACTUAL TO MAXIMUM NUMBER OF ARCS VS RE 


FORMAT 

FORMAT 

FORMAT 


('O'//13X, 
('O'//13X, 


FORMAT ('0'//13X, 

SIDUAL ERRORS') 

FORMAT C0'//13X, 'RATIO OF NODES TO ARCS VS RESIDUAL ERRORS') 
FORMAT C0«//13X, 'RATIO OF LOOPS TO ARCS VS RESIDUAL ERRORS') 


C0«//13X, 
(CONTINUED ON PAGE 3) 
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I 

FORMAT ('O'//13X,•EXECUTION TIME VS RATIO OF ACTUAL TO MAXIMUM NUM 
BER OF ARCS ' ) 

FORMAT ( 'O'//13X t * EXECUTION TIME VS RATIO OF NODES TO ARCS') 

FORMAT ('0'//13X,'EXECUTION TIME VS RATIO OF LOOPS TO ARCS') 

FORMAT (*0'//13X, 'NUMBER OF INPUTS VS PERCENT OF ARCS TESTED') 

FORMAT (»0'//13X, 'RATIO OF ACTUAL TO MAXIMUM NUMBER OF ARCS VS PE 
RCENT OF ARCS TESTED* ) 

FORMAT ('0'//13X, 'RATI G OF NODES TO ARCS VS PERCENT OF ARCS TESTE 

D' ) 

FORMAT ('0'//13X » 'RATIO OF LOOPS TO ARCS VS PERCENT OF ARCS TESTE 

D* ) 

FORMAT C0'//13X, 'NUMBER OF INPUTS VS EXECUTION TIME') 

FORMAT C0'//13X, 'NUMBER OF INPUTS VS PERCENT RESIDUAL ERRORS') 

FORMAT ('O'//13X,'EXECUTION TIME VS NUMBER OF INSTRUCTIONS') 

FORMAT Cl', 59X, 'DATA SUMMARY •//62X, 'PERCENT', 3X, 'PERCENT', 
2X, 'EXECUTION', 3X, 'REPAIR', 4X, 'ACTUAL', 4X, 'NODES', 5X, 

' LOOPS' / 4X, 'RUN', 35X, 'ERRORS', 3X, 'RESIDUAL', 2X, 'RESIDUAL', 
4X, 'ARCS', 6X, 'TIME', 6X, 'TIME', 5X, 'TO MAX', 6X, 'TO', 8X, 

'TO'/2X, 

•NUMBER', 5X, 'NODES', 4X, 'INPUTS', 4X, ' INSTR.', 4X, 'SEEDED', 

4X, 'ERRORS', 4X, 'ERRORS', 4X, 'TESTED', 4X, '(SEC)', 7X, 

'(HRS)', 4X, 'ARCS', 6X, 'ARCS',. 6X, 'ARCS'/) 

FORMAT ('O', 2X, 13, IX, **', 5X , 13, 7X, F5.0, 3(5X, F5.0), 

2(4X, F6.2), 2X, 2F10.5, 3X, F5.3, 2(5X, F5.3)) 

FORMAT ('O', 13X, 'THE LAST SEED USED WAS', Ill, '.') 

READ IN THE TYPE OF SIMULATION TO BE RUN (SIMNUM) 

***READ (5,100) SIMNUM 
READ IN NUMBER OF INPUTS TO BE USED (MINPUT) 

***READ (5,100) MINPUT 

READ IN THE MEAN NUMBER OF INSTRUCTIONS BETWEEN EACH NODE (MEANLN) 
***READ (5,100) MEANLN 

READ IN THE NUMBER OF GRAPHS TO BE EVALUATED (NUMOUT) 

***READ (5,100) NUMOUT 

READ IN WHETHER OR NOT THE NUMBER OF INSTRUCTIONS IS TO BE READ IN 
*.**READ (5,100) INREAD 

READ IN WHETHER DELETING AN ARC FROM THE STRUCTURE OR 
ADDING AN ARC TO THE STRUCTURE. 

***READ (5, 100) DELADD 

READ IN THE MEAN TIME TO EXECUTE AN INSTRUCTION IN MILLISECONDS 

I 

(CONTINUED ON PAGE 4) 
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294 
29 5 

296 

297 

29 8 

299 

300 

301 

30 2 
303 


304 
30 5 
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I 

(I TIME) 

***READ (5,100) ITIME 

READ IN MEAN TIME TO REPAIR AN ERROR IN MIN (MMTTR) 

***READ (5,100) MMTTR 

READ IN MEAN NUMBER OF INSTRUCTIONS THAT ARE ERROR FREE (MEANER) 
***READ (5,100) MEANER 

READ IN MEAN NUMBER OF ITERATIONS FOR EACH LOOP (MEANIT) 

***READ' (5,100) MEANIT 

READ IN MEAN NUMBER OF INSTRUCTIONS CHANGED WHEN AN ERROR IS 
FOUND (AVCHNG) 

***READ (5,100) AVCHNG 

READ IN THE ASTERICKS TO PUT AROUND OUTPUT ARRAYS (ASTER) 

***READ (5,206) ASTER 
READ IN THE NUMBER OF NODES IN) 

***READ (5,100) N 

READ IN THE ADJACENCY MATRIX (NODES) 


+ - + 

+ DO + 

+ 4- I = I , N + 

+ +-+ 


20 ++-m- + ***READ (5,140) (NOD ES ( I , J ) » J = 1,N) 


* 

* 


❖ IF * 

SIMNUM.EQ.l 



* 

& 



F 


(CONTINUED ON PAGE 5) 
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CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PAGE 


* 

* 


❖ 


* 

* 


% 


. * * . 

* IF * 

SIMNUM.EC.2 

1 

F 


. * * . 

# IF & 

SIMNUM.EQ.3 

* . . * 


IF ^ 

SIMNUM.EQ.4 


* 


* . jje 

! 

F 


. * * . 

. * IF * 

* SIMNUM.EQ.5 

* • 

* . * 

1 

F 


. * * . 

* IF * 

INREAD.EC.1 


* 


* . * 

I 

F 


* 

*- 


* 

*- 


* 

*■ 


❖ 

*- 


*- 


I WRITE(6 * 152) 


I WRITE(6,153) 


I WRITE(6,154) 


j WRITE!6,155) 


! WRITE(6,156) 


***WRITE (6,207) ITI ME, MMTTR, MINPUT, MEANER, AVCHNG 
***WRITE (6,208) MEANIT 
(CONTINUED ON PAGE 6) 
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I 

***WRITE (6,209) MEANLN 
***WRITE (6,204) NUMOUT 


PAGE 


NNQUT = 1 

IX =71286223 
IW = IX 


DETERMINE THE MAXIMUM NUMBER OF ARCS IN THE ADJACENCY MATRIX 

\_ _ 

I MAX ARC =N*(N-1) "* I 

NUMPTS IS THE LABEL FOR THE OUTPUT ARRAY 

I 

+■-+ 

+ DO + 

+ + «- + ++ + + + + 57 + 

+ + I = 1»N + 

+-- + 


| NUMPTS(I) =1 | 

- 

ASTER AND ASTOUT ARE TQ PUT ASTERICKS AROUND THE ARRAY 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ ---- 

+ I 

57 ++*+++++CQNTINUE 

CALL DIGRAF 

I 


I A$TOUT(I) =aster 


59 


■M- + + + + + + -M- 
+ 

•I- + 

+ 

+ 

+ + 

+ + 

++++++++++ 

+ + 

+ +■ 

+ 

+ 

+ 

* 

•K 

+ 

+ 

+ 


DO 

62 

I = 1,N 


1 


DO 

61 

J = 1,N 



iSEED IS THE ARRAY OF ERRORS SEEDED 

_ \_ 

J ISEED(I,JJ =0 


61 + f+ + + +CGNTINUE 

+ 

62 + + -M- + + + +CQNTINUE 


I 


(CONTINUED ON PAGE 7) 
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+ + + + + 


PAGE 


7 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

I 


SEED THE PROGRAM WITH ERRORS 


* 

* 


# IF # 

SIMNUM.EQ*I 

* . . * 

* . * 

i 

F 


*- 


| 1X^1722632 


CALL SEED 

THIS SEED CAN BE CHANGED TO THE LAST VALUE OF THE SEED ON THE 


63 


PREVIOUS TEST IF A CONTINUATION OF THE DATA IS DESIRED. 


* 

* 


. * * . 

* IF * . 

SIMNUM.EC.5 * 

. *- 

* . . * 

* . * 


| READ!5,108)IX 


ARCS =0, 
INST =0 
N S E E D = 0 


+ - 

+ DO 

+ + + + + + + + *-+ 65 

*■ + 1 = 1,N 

+ +- 

+ 

+ 

+ +-+ 

+ + DO + 

+ + -H-++-H- ++ 65 + 

+ + J=1,N + 

+ +- - -+■ 

+ I 

+ 

+ ZERO THE MATRIX WHICH RECORDS WHETHER OR NOT AN ARC HAS BEEN TEST 

+ 

+ I 


f TESTED( I » J) =0 


(CONTINUED ON PAGE 8) 
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8 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc.ee 

+ 

♦ 

+ 

- + 

+ . * * . 
t * ’nodes(iCj).neTi* * 

+ * . . *-7-7 

+ * . . * T | 65 | 

+ * , * - 

♦ I 

+ F 

+ 

+ 

+ 

+ CALCULATE THE NUMBER OF ARCS IN THE PROGRAM 

+ 

+ I 

+ ----- 

+ | ARCS = ARC S +1. I 

+ --- 

+ I 

+ CALCULATE THE NUMBER OF INSTRUCTIONS 

+ 

+ I 

+■ * — ” —-- — — “ “ * * — — —* 

+ I INST =INST+X(I * J) I 

+ ---.- 

+ I 

+ 

+ CALCULATE THE NUMBER OF ERRORS SEEDED 

♦ 

+ I 

+ — — — — — — — — ^ — — — ——— — — — — — — — — — — —-—— — — — 

+ 1 N$EED=NSEEO+ISEED{I,J) ! 

+ --- 

+ I 

4 - 

65 ++++++++CONTINUE 

CALCULATE THE NUMBER OF LOOPS IN THE PROGRAM 

_L 

I LOOPS=Q | 


+- + 

+ DO + 

++++++++++ 68 + 

+ + I=l,N + 

+ +- -f 

: i 

+ +—--+ 

+ + DO + 

+• + + + + *+4-* + 68 + 

+ + J=ltl + 

+ t-4- 

I 


(CONTINUED ON PAGE 


9) 
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CCCCCCCCCCCCCCCC.CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

+ 


• * ' * . 

$ IP % 

’nODES(I,J) • EQ* 1 


68 


4 - 
4 - 
4 - 
4 - 
+ 

4 - * 

+ * 

+ 

+ 

+ 

+ 

4 - 
4 “ 

4 - 

++++++++CONTINUE 

CHOOSE SAMPLE INPUT 


* 


I LOOPS--LOOPS 4-1 


I INPUT =1 


I 


77 


TOTREP =0• 

TIME =0. 

AVLINK =0. 

TTIME(NNOUT) =0, 
TREP(NNOUT) =0. 
TLINK(NNOUT) =0- 


78 


IZ = IX 
NODE =1 


79 


I 

***WRITE (6,227) 

NUMSUC IS THE NUMBER OF SUCCESSORS 

_l_ 

T NUMSUC =0 


+- 

4 - DO 

80 


+ 

+ 

+ 

+ 

+ 

4 - 

4 - 

4 - 

4 - 

4 - 

+ 

+ 

4 - 


+ J - I , N 
+- 


. * * . 

. * IF * . 

*NODES(NODE,JJ.EQ.l 

* . 

* . . * 

* . * 




* 

#- 


NUMSUC = NUMSUC 4-1 


80 4-4-4-4-4-4-4-4-CONTINUE 


(CONTINUED ON PAGE 10) 
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PAGE 10 


. * * . 

. * IF * . 

* NUMSUC.EQ.0 * 

* . . *- 

* | * 

F 


* 


. * * . 

# IF * 

NUMSUC.NE.l 


# . * 

! 

F 


* 

*- 


I 96 | 


82 


I K 


= 1 


I 83 ! 


82 


83 


THE SUCCESSORS ARE CHOSEN RANDOMLY WITH A UNIFORM DISTRIBUTION 
CALL RANDOM( IZ, U, I) 


\ 

K 

=1+NUMSUC*U 


1 

1 

r 

KK 

= 0 


1 


+ - 

+ DO 

* + + «• + ++++ 85 

+ J = I t N 
+--- 


(CONTINUED ON PAGE 11) 
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+ 
+ 
+ 
+ 
+ 
'+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ . 
+ 
+ 
+ 
+ 

+ 

+ 

+ 

+ 

* 

+ 

+ 

+ 


. * * . 

s(c J p £ ' 

*NODES! NODE, J KEQ.i * 
* . . *■ 

* . . * 

* . * 

I 

F 


* 

* 


IF 

KK.EQ.K 


I 

F 


❖ 


* 

#- 


| KK =KK+ 1 


86 


85 ++++**++CONTINUE 


86 


I L 


= J 


i 

TIME IS ASSUMED TO HAVE A MEAN OF IT IME MILLISECONDS PER 
INSTRUCTION WITH AN EXPONENTIAL DISTRIBUTION 


CALL EXPON(IZ, XTIME, 1) 
| XTIME=XTIME*ITIME 


* ’ * nodeTle.l * * * 

* . . *- 

* . . * 

* . * 


88 


CONVERT TIME INTO HOURS 

I 


TIME =TI ME +XTIME*ITER!NODE,L)*<X(NODE,L) 
+ X(L t NODE) )/ | 

I.0E03 


(CONTINUED ON PAGE 12) 
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PAGE 12 


I 89 | 


es 

89 


90 


_]_ 

T TIME =TIME+X(N3DE»L5*XTIME/1»0E03 I 

- 

***WRITE (6,230) N0DE,L 
TESTED(NODEiL) = 1 
IS THERE AN ERROR IN THIS PATH 

I 

. * * . 

. * IF * . 

*1 SEED(NODE,L i.EQ.O * 

* . . *- 

* . . * T | 95 I 

- 4 =.* - 

: 

ERROR FOUND 

WRIT E (6,2455 

SEE IF A NEW ERROR SHOULD BE INSERTED 


CALL NUSEED 


j ISEED(NODE,L) =ISEED(NODE,L5-1 

~r 


i 


CALL EX PON(IZ, U, 1) 

MEAN TIME TO REPAIR HAS AN EXPONENTIAL DISTRIBUTION WITH A 


MEAN OF MMTTR MINUTES 


| MTTR =MMTTR*U I 

„ .... 

TCTREP IS THE TOTAL REPAIR TIME IN HOURS 


I TOTREP =T0TREP+MTTR/60. | 


(CONTINUED ON PAGE 13) 
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I 73 f 


I 

*—• -r •— — — — — — — — — — — — — — — 

95 I NODE = L I 


I 79 ] 

— — 


AVLINK IS THE AVERAGE ARC T-RAVERSAL TIME IN MINUTES 


96 | AVLINK =(TIME/ARCS) I 

***WRITE (6,240) INPUT, TIME, TOTREP, AVLINK 

. * * . 

• * IF * . 

* SIMNUM.EQ.2 * 

Vi #____ _ _____ 

# " T | 9 7 | 

* . * -— 


TTIME(NNOUT) =TTI ME(NNOUT)+TIME 
TREP(NNOOT ) =TREP<NNOUT)+TOTREP 
TLINK(NNOUT) =TLINK(NNOUT)+AVLINK 


I 99 i 


I 


97 f NOUT =NN0UT-1 


(CONTINUED ON PAGE 14) 


PAGE 13 
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PAGE 14 


• * * • 

Jjt IP Jjt 

* * NiMOUT.GT.l * * 

* * * . . * T i 98 T 

# . * ---- 

I 

F 


TLINK(I) =AVLINK 
TREP(l) = T G T R E P 
TTIME(l) =TI ME 


\ 99 f 


98 I TLINK(NNCUT) =TLINK(NCUT) +-AVLINK 

TREP(NNOUT) =TREP(NOUT)+TOTREP 
1 TTIME(NNOUT) =TTI ME(NOUT)+TIME 


99 


TIME =0. 
TQTREP =0. 
AVLINK =0. 

INPUT =INPUT+1 


I 

KOUNT IS A COUNTER OF THE NUMBER OF SEEDED ERRORS REMAINING 


I KOUNT = O 


+ — 

+ DO 

++++++++++ 101 
+ + I = 1 T N 

+ + — 

+ 

+ 

+ + — 

+ + DO 

++++++++++ 101 
+ + J = 1,N 

+ +- 

I 




+ 

+ 

+ 

+ 


+ 

+ 

+ 


(CONTINUED ON PAGE 15) 


* 
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PAGE 15 


, * IF * 

NODES(I -J).EQ.l 

’ * . . * 

* | ^ 

F 


| KOUNT= KOUNT+ISEED(I ,J) 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

+ * NGDES< I.J) .EQ.i * 

+ * ♦ • *- 

+ 

+ 

+ 

101 ++++++++CONTINUE 

***WRITE (6,250) KGUNT, (NUMPTS(I)i I = 1»N) 

***WRITE (6,255 ) (A STOUT(I), I = i,N) 

1 

+-+ 

+ DO + 

++++++++++ 110 + 

+ + I = I t N * 

+ +-*-+ 

+ i 

110 +++++***WRITE (6,260) NUMPTS(I), (ISEED(I,J), J = 1,N) 

DISTINCT SEED FOR RANDOM NUMBER GENERATOR FOR EACH INPUT 

.L 

I IX =IX-12345 T 


* 

* 


* 

* 


* IF * 

SIMNUM. EC. 2 


* . * 

i 

F 


Ip 

SIMNUM.EQ.5*" 

* * 

F 


* 


I HI I 


* 

#- 


I ill ! 


MINPUT IS THE MAXIMUM NUMBER OF INPUT PATHS TO BE CHECKED 


(CONTINUED ON PAGE 16) 
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PAGE 16 


. * * . 

. * IF * . 

* INPUT.LE.MINPUT * 

* . . *- 

* . . * 

* . * 

I 

F 


I 78 | 


111 


CALCULATE THE NUMBER OF ARCS TESTED 


I TEST =0 


+- 

+- DO 

+ + + + + -M- + + 112 

+ I = IfN 
+- 


+ —- 

+ + DO 

++++++++++ 112 


+ J = 1,N 
+- 


<r 
+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

112 ++++++++CONTINUE 


. * * . 

* TESTED! I*»J) .EqTi* * 

* . . *- 

* . . * 


I TEST=TEST+1 


| INPUT=INPUT-l 
-- 

PTEST IS THE PERCENT OF ARCS TESTED 


| PTE ST(NNOUT) =100.*TEST/ARCS | 

~7 

F IS THE RATIO OF ACTUAL TO MAXIMUM NUMBER OF ARCS 

_L 

| F(NNOUT) = ARC S/MAXARC | 


(CONTINUED ON PAGE 17) 
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ccccccccccccccccccccccccccccccccccccccccccccceccccccccccccccccccc 

I 

FNA IS THE RATIO OF NODES TO ARCS 

I_ 

T FNA(NNOUT) =N/ARCS _ _ I 

- 

FLA IS THE RATIO OF LOOPS TO ARCS 

_L 

FLA(NNOUT) =LOUPS/ARCS 
TINST(NNOUT) =INST 
T INPUT ( NNOUT ) = INPUT' 

T SEED!NNOUT) =NS EE D 
TKOUNT(NNOUT) =KOUNT 
PSEED(NNOUT) =100.*(TKOUNT(NNOUT) 

/NSEED) 


PAGE 17 


I 

***WRITE (6,265) INPUT, INST, NSEED, KOUNT, PSEED(NNOUT), 
PTEST(NNOUT) 

***WRITE (6,270) TTIME(NNOUT), TRE P(NNOUT), T LINK(NNOUT) 
***WRITE (6,275) F(NNOUT), FNA(NNOUT), FLA(NNOUT) 

I NNOUT = NN0UT + 1 T 


* 

* 


* 

* 


. * * . 

► * IF * , 

NNOUT « GT•NUMOUT 


* , * 

i 

F 


* IF * 

SIMNUM.EQ.1 


* 

*- 


I 118 | 


* 

*- 


MEANLN=MEANLN+1 


(CONTINUED ON PAGE 18) 
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PAGE 




• 

* 

* . 



. * 

IF * . 


* 

SIMNUM.EC.l * 


* 

• 


. *■ 



* . 


. * 




* 

1 

I * 



• 

* 

* 

• 



. * 

IF * . 


* 

SIMNUM.EG.2 * 


* 

• 


. *■ 



* . 


• * 




* , 

. * 

1 




f 

L 



• 

* 

* • 



• * 

IF * . 



SIMNUM.EQ.2 * 


❖ 

• 


. *- 



* - 


• * 




* I 

r * 




F 

L 



+- 


-+ 



+ DO 


+ 

++++++++++ 113 

+ 

+ 


+ I 

« li N + 

+ 


+- 



+ 



1 

i 

+ 



1 


+ 


+- 



+ 


+ DO 


+ 



+ 113 

+ 



+ J 

= 1 

. ? N + 

+ 


+- 




113 


+ 

+ 

+ 


I ISEED(I,J|_=SVSEED(I , J) 
(CONTINUED ON P .AG E 191 


I 57 | 


1 INPUT=INPUT + l 


I 78 | 


l 


18 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc page 


115 


116 


+ 


* 

❖ 


* 


* 

* 


• * * . 

$ J P 3{( 

SIMNUM.EQ.5 

« . . * 

* . sf= 

I 

F 


• ^ 

* IF * 

DELADD.EQ.O 

* . . * 

I 

F 


. * * . ■ 

* IF * 

DELADD.EQ.1 

« . . * 

* J * 

F 


* 


* 

*- 


* 

*« 


63 


1 CALLDELARC 


I CALLADDARC 


***WRITE (6,180) (NUMPTS(I), 1=1,N) 
***WRITE (6,255) (ASTOUT(I), I = 1,N) 


+ -- 

+ DO 
117 
+ I = 1, N 

+ --- 

I 


117 +*+++***WRITE (6,260) NUMPTS(I), (NODES(I,J), J=1,N) 
***WRITE (6,203) 


i 63 f 


118 ***WRITE (6,305) IZ 

***WRITE (6,284) 

CALL PLOTP (TTIME, 

***WRITE (6,285) 

***WRITE (6,204) 

I 

(CONTINUED ON PACE 20) 


TKOUNT, NUMOUT, 0) 


19 
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CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PAGE 20 

I 

CALL PLOTP (TTIME, PSEED, NUMOUT, 0) 

***WRITE (6,286) 

***WRITE (6,234) 

CALL PLOTP (TREP, TKOUNT, NUMOUT, 0) 

***WRITE (6,287) 

. * * . 

. * IF * . 

* SIMNUM.NE.2 * 

3jC _____ 

* . . * T | 119 | 

* | * - 

F 

I 

***WRITE (6,284) 

CALL PLOTP (TINPUT, TKOUNT, NUMOUT, 0) 

***WRITE (6,288) 

***WRITE (6,284) 

CALL PLOTP (TINPUT, PTEST, NUMOUT, 0) 

***WRITE (6,296) 

***WRITE (6,284) 

CALL PLOTP (TINPUT, TTIME, NUMOUT, 0) 

***WRITE (6,300) 

***WRITE (6,284) 

CALL PLOTP (TINPUT, PSEED, NUMOUT, 0) 

***WRITE (6,301) 

119 . * ^ . 

. * IF * . 

*( SIMNUM.EQ.2).OR.(SIMNUM.EQ ,5 ) 

# o £___ 

* . . * T | 121 | 

J|S g „_ 

F 

I 

***WRITE (6,284) 

CALL PLOTP (TTIME, TINST, NUMOUT, 0) 

***WRITE (6,302) 

***WRITE (6,234) 

CALL PLOTP (TINST, TKOUNT, NUMOUT, 0) 

***WRITE (6,269) 

(CONTINUED ON pIgE 21) 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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. * * . 

. ♦ IF * . 

* SIMNUM.EQ.l ♦ 

♦ . • *- 

* . . * T I 121 I 

* . * - 

I 

F 


♦♦♦WRITE (6,234) 
CALL PL OTP (F, 
♦♦♦WRITE (6,290) 
♦♦♦WRITE (6,284) 
CALL PLOT P (FN A, 


♦♦♦WRI TE t 6,291 ) 
♦♦♦WRITE (6,284) 

CALL PLOTP (FLA, 
♦♦♦WRITE (6,292) 
♦♦♦WRITE (6,284) 

CALL PLOTP (TTIME, 
♦♦♦WRITE (6,293) 
♦♦♦WRITE (6,284) 

CALL PLOTP (TTIMEt 
♦♦♦WRITE (6,294) 
♦♦♦WRITE (6,284) 

CALL PLOTP (TTIME, 
♦♦♦WRITE (6,295) 
♦♦♦WRITE (6,234) 

CALL PLOTP (F, 
♦♦♦WRITE (6,297) 


♦♦♦WRITE (6,284) 
CALL PLOTP (FNA, 
♦♦♦WRITE (6,298) 
♦♦♦WRITE (6,234) 
CALL PLOTP (FLA, 
♦♦♦WRITE (6,299) 


(CONTINUED ON PAGE 22) 


TKOUNT, NUMOUT, 0) 


TKOUNT, NUMOUT, 0) 


TKOUNT, NUMOUT, 0) 


F, NUMOUT, 0) 


FNA, NUMOUT, 0) 


FLA, NUMOUT, 0) 


PTEST, NUMOUT, 0) 


PTEST, NUMOUT, 0) 


PTEST, NUMOUT, 0) 
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121 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


+- -+ 

+ 00 + 

125 + 

+ K = 1, 2 + 

+-'- 

I 

***WRITE (6,303) 

***WRITE (6,255) (ASTOUT(I), I = l,N) 

+ ----+ 

+ DO + 

122 + 

+ I = 1,NUMOUT + 

+-+ 


122 +++++***WRITE (6,304) I, N, TINPUT(I), TINST(I), TSEED(I), TKOUNT(I), 
+ PSEED(I), PTEST(I), TTIME ( I ) , TREP(I), F(I), FNA(I), FLAU) 

125 ++++++++CONTINUE 


. * * . 

. * IF .* . 

♦(SIMNUM.NE.3).AND.< S IMNUM.NE.4) 

* . , 4 - 

* . . * T I CALLGRAPHO 

* . # 



NEED A MINIMUM OF 10 GRAPHS TO ANALYZE BEFORE CALLING HISTG 

I 

. * * . 

. * IF * . 

* SIMNUM.NE.5 . * 

* . . * T I 126 | 

* $ _ _ _ 

F 

I 

CALL HISTG (TTIME, NUMOUT, 0) 

***WRITE (6,280) 

CALL HISTG (TREP, NUMOUT, 0) 

***WRITE (6,281) 

CALL HISTG (PSEED, NUMOUT, 0) 

***WRITE (6,282) 

CALL HISTG (PTEST, NUMOUT, 0) 

***WRITE (6,283) 

126 CONTINUE 

, I 

(CONTINUED ON PAGE 23) 
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PAGE 


1 


100 
130 
140 
150 
160 
* 130 

190 
V 200 

210 

222 

225 

230 


SUBROUTINE DIGRAF 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

CONSTRUCT THE DIRECTED GRAPH C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COMMON/ALLL/N, N0DES(30,30) 

COMMON/ARRAY/X( 30 » 30 ) , NUMPTSC30), AST0UT(30), AVC.HNG 
COMMON/GEN/ME ANLN, IX 

CGMMON/ERROR/ISEED(30,3O) , ITER<30,30), MEANER, INST, MEANIT, IZ 

COMMON/GRAF/INREAD 
INTEGERS CHANGE, AVCHNG 
FORMAT (16F5.0) 

FORMAT («0*, 10X, 14, IX, 'NODES IN DIRECTED GRAPH.'//) 

FORMAT(* 0 *, 13X, 'MATRIX OF ARC LENGTHS‘//7X , 30(14)) 

FORMAT ('O', 15, • *', 30(F4.0)i 

FORMAT ('O', 13X, 'NOTE: 0.0 INDICATES NO ARC FROM I TO J'//) 
FORMAT ('O', 13X, 'ADJACENCY MATRIX'//7X, 30(14)) 

FORMAT ('O', 15, ' *'» 30(14)) 

FORMAT ('O', 13X, 'NOTE: NONEXISTENCE OF AN ARC IS INDICATED *3Y A 
ZERO'//) 

FORMAT (• ', 7X, 30(A4)) 

FORMAT (• », 7X, 30(A4)) 

FORMAT ('O', 15, ' **, 30(14)) 

FORMAT ('O', 13X, 'LOOP ITERAT IONS*//7X, 30(14)) 

IF INREAD = 0 RANDOMLY GENERATE INSTRUCTION LENGTHS 

IF INREAD = 1 READ IN INSTRUCTION LENGTHS 
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PAGE 




+- - + 

+ 00 + 

10 + 

+ + I = 1,N + 

+ +---+ 

+ I 

10 +++++***READ (5,100) (X(I,J), J 

T 25 T 


= 1,N) 


GENERATE THE NUMBER OF INSTRUCT IONS' BETWEEN EACH NODE - 


EXPONENTIALLY DISTRIBUTED 


20 


+ - 

+ DO 

+ + + {■ ^+ + + f + 2 5 


—+ 


+ 

+ 

+ 

+ 

+ 

4 - 

+ 

4* 

+ 

+ 

4 - 

4 * 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4 - 

+ 

4 - 

+ 

+ 

4 - 

+ 

■f* 

+ 

+ 

+ 

+• 

+ 

+ 

4 - 


+ 1 = 1,N 

+- 


+- 

+ DO 

25 

+ J=1,N 

4 -- 


| X(I,J1 =0.0 


* 

* 


* IF * ■ 

NGDESd, J» • EQ. 0 


* 


I 25 | 


CALL EXPON ( IX,XLNTH, 1) 

MEANLN IS THE EXPONENTIALLY DISTRIBUTED MEAN LENGTH OF EACH ARC 

|____ 

I X(ItJ) =MEANLN*XLNTH+1.0 I 

T 


25 ++++++++CONTINUE 

***WRITE (6,130) N 
***WRITE (6,140) (NUMPTS(I), I 
***WRITE (6,210) (ASTOUT(I), I 
(CONTINUED ON pIgE 3) 


= 1, N) 
= 1,N) 
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PAGE 3 


+-* 

+ 00 + 

++++++++++ 30 + 

+ + 1=1,N + 

+ + —-+ 

+ i 

30 44*44***WRITE (6 v 150 j NUMPT 511), (X(I,J), J=1,N) 
***WRITE (6,160) 

***WRITE (6,180) (NUMPTS(l), 1=1,N) 

***WRITE (6,210) (AST OUT(I), I 

- -+ 


It N) 


4 

4 

4 

4 


4 - 

4 DO 

50 

4 1=1,N 

4- 


I 


50 44*44***WRITE (6,190) NUMPTS(I), (NODES(I,J), J=1,N) 
***WRITE (6,200) 

THE NUMBER OF ITERATIONS PER LOOP IS DISTRIBUTED WITH / 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


DISTRIBUTION WITH A MEAN OF MEANIT 

I 

4--4 

4 DO 4 

55 4 

4 I = 1,N 4 

4-4 


4--- 

4 DO 

55 

+ J = 1 , N 

4- 


i iT er ( i » j J =o 1 

. * 

* . 

. * IF * . 

NODES( I . 

r J) • EQ. 0 * 

• 

• *----*- 

* . 

. * T | 55 | 

* , 

1 

. * - 

1 ITER(IyJ) =L | 


(CONTINUED ON PAGE 4) 


UNIFORM 
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PAGE 


4 


. * IF * 

* J.GT.I 

* . 

# . . * 

* . * 

I 

F 


CALL RANDOM (IX, TITER, 1) 


I I TER(I , J) =1+MEANIT*TITER 


55 ++++++++CONTINUE 


55 | 


CREATE LABELS FOR THE OUTPUT ARRAY 


***WRITE (6,230) (NUMPTS(I), 1=1,N) 
***WRITE (6,222) (ASTGUT(I), 1=1,N) 

+ + D0 *+ 

H-+ + + + 73 + 

+ I = 1, N + 

+ - + 


73 +++++***WRITE (6,225) NUMPTS(I), (ITER(I,J), J=1,N) 
RETURN 
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100 

110 


SUBROUTINE ADDARC 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

ADDARC ADDS AN ARC BETWEEN TWO NODES C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COMMON/ALLL/N, NODES(30,30) 

COMMON/ARRAY/XC30 ,30}, NUMPTS(30), AST OUT(30) , AVCHNG 
C GMMON/GEN/ME ANLN, IX 
FORMAT (215) 

FORMAT (*1», 13X , 'A NEW ARC HAS BEEN ADDED FROM NODE', 13, 

• TO NODE', 13, • .OF LENGTH', F4.0) 

ADD A NEW ARC FROM N1 TO N2 
***READ (5,100) Nl, N2 

| NODES(Nl,N2) =1 T 

- 

CALL EXPON (IX, XL NTH, l) 

MEANLN IS THE MEAN LENGTH OF THE ARC TO BE EXPONENTIALLY DISTRIB. 

__ 2 _ __ __ ____ 

I X(N1,N 2) =MEANLN*XLNTH+1.0 i 


***WRITE (6,110) Nl, N2, X(N1,N2) 
RETURN 









PAGE 


1 


SUBROUTINE OELARC 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

DELARC DELETES AN ARC BETWEEN TWO NODES C 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMMON/ALLL/N, N0DES(30,30) 

COMMON/ARRAY/X(30,30), NUMPTS(30), AST OUT( 30) , AVCHNG 
COMM0N/ERROR/ISEED<30,30), ITER(30,30>, MEANER, INST, MEANIT, IZ 

CCMMON/OUT/NSEED, INPUT, SVSEED(30,30) 

INTEGER SVSEED 
100 FORMAT (215) 

110 FORMAT Cl*, 13X, 'AN ARC WAS DELETED FROM NODE', 13, 

• TO NODE* , 13, ' . 

***REAO (5,100) Nl, N2 

N0DES(N1,N2) =0 
X(N1» N2 ) =0.0 
ISEED(N1^N2) =0 

T 

SVSEED(N1, N 2) = 0. 

*#*WRITE (6,110) Nl, N2 
RETURN 


END 
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PAGE 


1 


♦ 

1 

\ 


■* 


* 

SUBROUTINE SEED 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCGCCCCC cc 

c 

SEED THE GRAPH WITH ERRORS C 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccxccccccccccccc 

205 

COMMON/ALLL/N, NODES<30,30) 

COMMON/ARRAY/X(30,30) , NUMPTSI30), AST0UT(30), AVCHNG 

COMMQN/ERR0R/ISEED(30,30), ITER(30,30), MEANER, INST, MEANIT, IZ 

CCMMQN/OUT/NSEED, INPUT, SVSEED ( 30,30 ) 

COMMON/GEN/MEANLN, IX 

INT EGER*4 CHANGE, AVCHNG, SVSEED 

DIMENSION ER(30), IER(30) 

FORMAT ('O', 13X, 'SEED ERRORS AT INSTRUCTION INTERVALS'// 10X, 
3014) 

210 

FORMAT Cl'// 13X, 'TOTAL NUMBER OF INSTRUCTIONS IS', 15) 

22 0 

FORMAT CO'// 13X, 'TOTAL NUMBER OF ERRORS SEEDED IS', 14// 13X, 
•THE ERROR MATRIX'//7X, 30(14)) 

22 2 

FORMAT (• ', 7X, 30(A4)) 

225 

FORMAT ('O', 15, • *• , 30(14)) 

GENERATE ERRORS WITH EXPONENTIAL DISTRIBUTION 

CALL EXP ON(IX, ER, N) 

SPREAD ERRORS TO GIVE MEAN OF MEANER 


1 

+- + 

+ DO + 


+*++++++++ 60 + 


+ 

t 

* + 
¥ 

¥ 

¥ 

¥ 

¥ 

¥ 

+ I = 1,N + 

+*-- ¥ 

1 

1 IER(I) =MEANER*ER<I) 1 

| 92 

(CONTINUED ON PAGE 2) 








PAGE 


2 


60 + + + + + + +■+CQNTINUE 


in =1 


NSEED IS THE NUMBER OF ERRORS SEEDED 


| NSEED=0 I 

____ 

INST IS THE NUMBER OF INSTRUCTIONS 


| INST =0 


NUMER IS THE ERROR SEED/INSTRUCT I ON COMPARATOR 


I NUMER=0 I 
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PAGE 


3 


+ 

+ 

+ 

4 - 

+ 

+ 

+ 

+ 

4 - 

+ 

♦ 

4 - 

4 * 

4 - 

+ 

+ 

+ 

4 * ■ 

63 + 


4 * 

4 - 

4 - 

+ 

4 - 

4 - 

4 * 

+ 

4 - 

+ 

4 - 

4 - 

+ 

+ 

+ 

4 - 

+ 

4 - 

4 - 

+ 

4 - 

4 * 

+ 

+ 

4 - 

+ 

+ 

+ 

+ 

4 - 

+ 

4 - 

4 - 

+ 

+ 

+ 

4 - 

+ 

4 * 

4 - 

4 * 

4 

4 

4 

4 - 

4 

64 4 


THROUGH THE LOOP 


. * * 
* IF 
J.LE.I 


# . # 

i 

F 


I IER< II ) = IER< I Ii/( < X < I , J)*ITER< I , JM 


NUMEK=NUMER+IER<11) 


. * * . 

Ip Jjt 

INST.LT.NUMER 


I SEED <I»J) =ISEED(I»J)+l 
NSEED=NSE ED + 1 

II =11+1 


. * * . 

* IF * . 

II.LE.N * 

. *- 

* ^ * 

F 


I 63 | 


II =1 


I 63 | 


| NUMER=NUMER-IER< 11) 


(CONTINUED ON PAGE 
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PAGE 


4 


6 5 + <■ + + + ■*■ + + CON TINU E 
+ 

70 ++++++++C0NTINUE 


+ - 

+ DO 

71 

+ I = i * N 
+- 


+ -- 

+ DO 

71 

♦ J = 1*N 


+ SVSEED IS THE SAVE COPY OF THE ORIGINAL SEEDED ERRORS 

+ 


71 + + + + + + + 4-I SVSEED(IjJ) =ISEED(I,JI 


NSEED IS THE NUMBER OF ERRORS SEEDED 


***WRITE (6,210) INST 

***WRITE (6,205) (IER(I), I = 1,N) 

***WRITE (6,220) NSEED, (NUMPTS(I), I = 1,N) 
***WRITE (6,222) (ASTOUT(I), 1=1,N) 


+- 

+ DO 

75 

+ I = 1, N 
+- 


***WRITE (6,225) NUMPTS(I), (ISEED(I,J), J=1,N) 


75 ++++++++C0NT1NUE 


RETURN 
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PAGE i 


SUBROUTINE NUSEED 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 


CORRECTING PREVIOUS ERROR MAY CREATE A NEW ERROR 


C 


c 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMMON/ALLL/N, N0DES(30,30) 

COMMQN/ARRAY/X(30♦30)» NUMPTS(30), AST0UT<30), AVCHNG 

COMMON/E RROR/ISEED(30»30), ITER<30,30), MEANER, INST, MEANIT, IZ 

CCMMON/OUT/NSEEO, INPUT, SVSEED(30,30) 

COMMON/NEW/IW 
INTEGERS CHANGE, AVCHNG 
DIMENSION E R R (2 ) 

FORMAT ('O', I3X, 'NEW ERROR INSERTED FROM NODE', 14,' TO NODE', 
14//) 


| NUM =0 


THE ERROR IS EQUALLY LIKELY TO RESULT IN CHANGES TO ALL ARCS 


CALL RANDOM (IW, ERR, 2) 



=1+N*FRR(1) 
= I+N* ERR(2) 


NEW.ERROR CREATED AT (I,J) 


IS THIS A NODE OF THE DIRECTED GRAPH 


. * . 

* ‘nodes( I , J) .EQ.I 


T | 20 


* , * 



= J + 1 

F(J.GT ,N)J =1 


(CONTINUED ON PAGE 2) 


9 6 














PAGE 


2 


20 


I 

NUM=NUM+ 1 

IF NODE I IS A TERMINAL NODE NO NEW ERROR IS INSERTED 


* 

* 



* IF * . 

NUM.GT.N * 

• *- 

* . . * T | RETURN 

* | * - 

F 


I 10 | 


CHANCE OF INSERTING A NEW ERROR IS BASED ON THE NUMBER OF 
INSTRUCTIONS THAT MUST BE CHANGED 
CALL RANDOM (IW, ERNU, 1J 

AVCHNG IS THE ASSUMED (INPUT) AVERAGE NUMBER OF INSTRUCTIONS 
CHANGED WHEN AN ERROR IS FOUND 

CHANGE IS THE UNIFORMLY DISTRIBUTED NUMBER OF INSTRUCTIONS CHANGED 

_L 

I CHANGE =1+AVCHNG*ERNU*2. T 

- 

CRATIO IS THE RATIO OF INSTRUCTIONS CHANGED TO THE NUMBER OF 
INSTRUCTIONS IN THE UNIFORMLY RANDOMLY CHOSEN ARC 
_ [_ 

| CRATIO =CH ANGE/Xl I ? J ) | 

_ 

THE CRITERIA FOR INSERTING A NEW ERROR IS UNIFORMLY DISTRIBUTED 
CALL RANDOM (IW, RATIO, 1) 

IF THE CHANGE RATIO (CRATIO) IS GREATER THAN THE CRITERIA (RATIO), 
THEN A NEW ERROR IS INSERTED 
(CONTINUED ON PAGE 3) 


97 















PAGE 


. * IF * . 

CRATIO.LE.RATIO * 

. . * --- 

* . . * T | RETURN 

* . * 

F 


I SEED( I.J) = 1 SEED{ I,J) + 1 
NSEE D=NSE E D+I 


***WRITE (6,247) I,J 


RETURN 













+ + 


PAGE 1 


I 


SUBROUTINE GRAPHO 

cccccxcccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


GRAPHICALLY DISPLAY THE DIRECTED GRAPH USING THE CALCOMP PLOTTER C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COMMON/A LLL/N, NODES(30,30) 

DIMENSION ND(30,30) f X(30), XX<2), Y(30), YY(2) 

DIMENSION XL(36), YL(36), THETA(36) 

INTEGER-4 I TB ( 12) / 12*0/ 

REAL-4 RTB<28)/28*0.0/ 

ITB IS DRAWP ARRAY 


4- 

+ 

4- 

+ 

4H 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 


ITB(l) =1 

ITB(2 ) =0 

ITB(3) =9 

ITB(4) =15 

ITB(8 ) =2 

I TB(9 ) =0 

ITB(IO) =2 

1TB(11) =0 

1 


IS DRAWP ARRAY 


1 


j RTB<1) =1.0 

7 

| RTB(2) =1.0 

1 

1 

4- 


+ DO 

4- 

+ + 4 

4- 

+ I = 1»N 

4- 

4- - 

| 

-4- 

1 

+ _- 

-4- 

+ DO 

4- 

+ + 4 

4- 

+ J = 11N 

4- 

+ -*- 

-4- 

1 


S ARRAY USED TO 

CHECK FOR PREDECESSORS 

1 


| ND(I,J) =0 1 


(CONTINUED ON PAGE 2) 
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PAGE 


■f 

4 + *++++++CONTINUE 


X(1) =4.5 

Y<1) =15.0 
X ( 2) =X(1) 

— «(«.«.«» — — — — — ■ —. _ _ — _ _ — ■ 

I 

CONST IS Y-DI RE CTION SCALE FACTOR 


C0NST=30./N 
Y(2J = Y(1)-CON ST 


I 

CALL DRAWP(2,X,Y,ITB,RTB) 


| XCONST =1.5 


KOUNT IS X-DIRECTION SHRINK FACTOR - SHRINKS AFTER EVERY OTHER 
FANOUT 

I 

I K0UNT=1 | 

-- 

LEVEL IS THE LEVEL IN THE DIRECTED GRAPH 


LEVE L= 2 
NN =N-1 
ITB(l) =2 


I 

LABEL IS THE NUMBER OF EACH NODE 


.1 LAB EL= 2 7 


+ - + 

+ DO + 

+++++++♦++ 5 + 

+ + M = 11N + 

+ +---+ 


(CONTINUED ON PAGE 3) 



















PAGE 


3 


* 

* 


KK IS THE LAST NODE IN THE GIVEN LEVEL 


+ 

4* 

4* 

4-- 
4* 

4* 

4- 
4- 
4- 
4- 
4- 
4- 
4- 
4* 

4- 
4- 
4- 
4- 
4- 
4-' 

4- 
4- 
4- 
+ 

4- 
4- 
4- 
4- 
4- 

4* * . 

4- * 

4- 

4- 

4- 

4- 

4- 

5 4-4- + + 4-4- + +CONTINUE 


# * # 

. * IF * ■ 
NODES(1,M) • EQ« 0 


& 


I 5 


KK 

KT 


= M 

= KK + I 


. * IF * . 

* NODESl1,KT).EQ.O 

* „ 


I 

F 


* 

*- 


I 6 


KM IS THE LAST NODE IN THE NEXT LEVEL 

I 

+■-4* 

4* DO 4- 


4- 

4- 

4- 

4- 

4- 

+ 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4 

+ 

4- 

4- 

4- 

4- 

+ 


4- M = KK f N 
4-- 


. * * . 

. * IF * . 

* NODES(KK j M).EQ.O * 

* • . 

* • • * 

* 


I 7 



(CONTINUED ON PAGE 4) 
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. * IF * . 

* NODES(KK,KX).EQ.O * 

* . . 

* . . * 

* * 

F 


;ONTINUE 

I 

+-+ 

♦ DO + 

- + 50 + 

+ I = 2 f NN + 

+-+ 


II =1-1 


NUMSUC IS THE NUMBER OF SUCCESSORS 


i NUMSUC =0 


+ --- + 

+ 00 + 

20 

+ J = I,KM + 

+---+ 


. * IF * . 

* NODES! I, Ji .EQ.O * 

* • . 

❖ . . *■ 

« „ * 


4 --+- 

+ DO + 

•++++++ 18 + 

+ K = 1, I I + 

+- - -+ 

I 

IS THERE A NODE ALREADY IN EXISTENCE 


(CONTINUED ON PAGE 5) 
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17 


* 

* 


• * 

* IF 
NO(K * J}, 


I 

F 


E G * 0 


* 


18 


+ 

+ 

+ 
f 
+ 

+ 

4- 
4- 
+ 

4- 
+ 

4- 
4- 
+ 

4- 
4- 
+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4- 
4- 
4- 
4- 
+ 

+ 

+ 

4- 
4- 
+ 

4- - 

4- 

4- 

4- 

4- 

+ + + 4-++ + +CONTINUE 
♦ I 

+ - 

+ | NUMSUC =NUMSUC +1 


* 


18 


CONNECT NEW NODE TERMINATING AT EXISTING NODE 

__ [ _ 

XX( 1 J =X( I ) 

XX(2)=X(J) 

V Y ( 1) = V ( I ) 

YY(2)=Y(J) 


1 

CALL DRAWP(2,XX t YY,ITB,RTB) 


I ND(ItJJ =1 


I 20 | 


20 


4- 
+ 

+ 

++4- + + + ++C0NTINUE 
+ 

+ 

+ 

4- 

4- * 

4- * 

4- 
4- 
+ 

4- 
4- 
4- 
4- 


. * * . 

j{c jp ^ 

NUMSUC.L E•0 


* 


! 44 | 


I 

F 


(CONTINUED ON PAGE 
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25 


+ 

+ 

4- 

4* 

4- 

4- 

+ 

4- 

+ 

4* 

+ 

4- 

4- 

4- 

+ 

+ 

4- 

4- 

+ 

4- 

4- 

4- 

+ 

+ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4* 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


* 

* 


♦ * * 

* IF * 

NUMSUC.NE.1 

* • . * 

* • * 

1 

F 


* 

*- 


I 30 | 


FANOUT TO ONE NEW NODE 

I 


L =LAB EL + 1 
LAB EL = LAB EL + l 
NO( I f L ) =1 
X(L) =X(II 

Y ( L) =Y tI )-CONST 
XXC1)=X <I ) 

X X ( 2 ) = X ( L ) 

YY( 1)=Y(I) 

Y Y { 2 ) = Y { L ) 


I 

CALL DRAWP(2tXXtYYtITBtRTB) 

IS THERE A NODE ABOVE THAT NEEDS CONNECTING TO THIS NODE 

I 

4--4 

4* DO 4- 

28 4- 

4 K = 1,11 4 

4 -~- —4 


❖ 

* 


* IF * 

ND(K * L)•E C • 1 

* a • * 

* * 

I 

F 


• * * • 

* * IF * . 

* NODES(K,L)• EQ. 0 * 

* # - 

* . • * 

i 

F 


28 


I 28 | 


XX( 1 ) = X lK ) 
Y Y ( 1 ) = Y ( K ) 


(CONTINUED ON PAGE 
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30 




+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


XX{2)=X < L 


! 


1 YY(2)=Y(L> 


T 

| ND(KtL) =1 


l 


i 

CALL DRAWP12 ,XX,YY,ITB,RTB) 


28 ++++++++CGNTINUE 

+ [ 

+ I 44 j 

+ --- 

+ 

+ 

+ 


+ 

+ 

+ 

+ 

+ 

* 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

♦ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+• 

+ 

* 

+ 


FANOUT TO MULTIPLE NODES 

I 

— — — — ^ — 

I FAN =2.*XC0NST/<NUMSUC-1 ) I 

-— — — — — — — —- — — i 

I 

+-+ 

+ DO + 

35 + 

+ L = 1» NUMSUC + 

+---+ 

I 

M = LAB EL + L 

ND(I,M ) =1 


* 

* 


. * * 

* IF 
L.GT.I 

* . 

I 

F 


* . 

* 

. *- 

* 


I 31 


| X(M) = X(I)-XCONS T 


I 33 j 


LM =M-1 
X(M) = X t LM)* FAN 


Y (M ) =Y(I)-CONST 
XX(1)=X(I) 


(CONTINUED ON PAGE 8J 
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YYl 1) 


=Y ( 11 


XX(2)=X(M) 
YY(2)=Y(M) 


CALL DRAWP(2,XX»YY,ITB,RTB) 

IS THERE A NODE ABOVE THAT NEEDS CONNECTING TO THIS NODE 


+ - 

+ DO 

34 

+ K = 1 , I I 
+- 


# IF ; 

ND(K,M).EQ.I 


* NODES(K,M) .EQ.O * 


! 34 


I 34 I 


XX(i)=X(K> 
YY(1J=Y(K) 
XX(2)=X(M) 
Y Y ( 2 J = Y { M ) 
ND(K,M) =1 


+ CALL DRAWPI2,XX,YY,1TB,RTB) 

+ 

34 + f- + n-+ + *CONTINUE 

35 +++++♦++CONTINUE 

+ I 

+ --- 

+ I LA6EL=LABEL+NUMSUC 

+ I K OUNT= KOUNT+I 

+ --- 


(CONTINUED ON PAGE 
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PAGE 
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40 


44 


+ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 ” 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


* 

* 


• * * 

* IF * 

KOUNT.LT.2 

* „ * 

i 

F 


* 

*- 


1 44 | 


REDUCE THE X-DIRECTIGN SCALE FACTOR 


XCONST =XCCNST/1.5 
K 0 U N T = 0 


IS THIS NODE STILL ON THE GIVEN LEVEL 


* 

* 


. * * . 

* IF * 

I .LT.KK 

* . * 

I 

F 


* 


50 


NO, CHANGE LEVELS 

I 


j LEVEL= LEVEL + I 

T i 

DETERMINE THE LAST NODE IN THE NEW LEVEL 


+- 

+ DO 

45 

+ M = I, N 

*■ -- 


• * * • 

t # j p # 

* 'nodes(I,M) . EQ.0 * * 

* . . *- 

* * 

F 


I 45 


I KK =M 


(CONTINUED ON PAGE 10) 
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KT = KK+1 


. * * 

# IF 

* NODE S(1,KT) . 

* . 

* . 

* . * 


45 


46 


I 


4 - 

4 DO 

47 

4 M = KK,N 

4 —- 


47 

50 


4* 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

44444444C0NTINUE 

4 DETERMINE THE LAST 

4 
4 
4 
4 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 * . 

4 * 

4 
4 
4 
4 
4 

44444444C0NTINUE 
4 

44444444C0NTINUE 

GRAPH LOOPS 


* < 

EQ.O 


* 

❖- 


I 46 


NODE IN THE NEXT LEVEL AFTER THE PRESENT LEVEL 


# * * 

. * IF * , 

* NODE S{KK * M).EQ.O 

* • 

* • . * 

* . * 


* 

*- 


KM 

KX 


= M 

= KM4 1 


# * * 

# j p 

* NODES( KK t KX) 

* . 


* . 

EQ.O * 

. *- 

* 


1 47 | 


I 50 | 


(CONTINUED ON PAGE 11) 
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I K 


* 

* 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ H 

++++++++ + + 
+ i 

+ 

+ 

+ 

+ 

+■ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

♦ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


DO 

90 

I = 2, 

»N 



= I- 

-1 



DO 

80 

J = 1, 

»K 


. * * . 

» * IF * « 

NODES { I ,J) • EQ. 0 


❖ 

*- 


80 


I 

F 


PI 


=3.1415926 


ESTABLISH THE DIRECTION OF CURVE FOR LOOP 


* 

* 


. - * . 

# i p 

X(I )•LT.X(J) 


I PI=-PI 


STARTING POINT FOP CURVE 

f XL(n=xm 

|^VL ( 1)= y ( i ) 

ENDING POINT FOR THE CURVE 


XL < 36) =X(J) 
YL(36) =Y(J) 


(CONTINUED ON PAGE 12) 
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XDIF AND YD IF ARE THE DIFFERENCE BETWEEN THE X AND Y VALUES RESP, 


XDIF =X(I )-X(J) 
YDIF =Y(I )-Y(J) 


R IS THE RADIUS OF THE LOOP 


j R =0.5*SQRT(XDIF**2+YDIF**2) 


< XC »YC) IS THE CENTER OF THE LOOP 


XC =(X(I)+X<J))/2. 
YC =(Y(I)+Y(J))/2. 


THETEdl IS THE INITIAL ANGLE 


I THETA(l) =ATAN2(YDIF,XDIF) 


CREATE 36 POINTS IN A SEMICIRCLE FORMING THE LOOP 


+- 

+ DO 

70 

+ L = 2,35 

+ - 


M =L-1 

THETA(L) =THETA(M)+PI/36. 
XL(L)=R*COS(THETA(L))+XC 
YLIL)=R*SIN(THETA(L))+YC 


70 ++++++++CONTINUE 
+ 

+ CALL DRAWPI36,XL,YL,ITB,RTB) 

+ 

80 ++*-++ + +^CONTINUE 

+ 

90 + + + +++ + +CONTINUE 

DRAW DIAMONDS AT EACH NODE 


ITB(l) =3 
ITB(2 ) =4 


CALL DRAWP(N ,X,Y,1TB,RTB) 
RETURN 


END 
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ERROR 5IMUIATIQW MODEL COMPUTER PROGRAM CODE 


w 


//GREEN05 JOB (1723,0409,CS42),'GREEN SMC 2827 *,TIME=4,MSGCLASS=0 
// EXEC FORTCLGP, REGION. G0=200K 
//FCRT.SYSPRINT DD SYS0UT=0 
//FORT.SYSIN DD * 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C SIMULATION OF ERRORS IN SOFTWARE 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DIRECTIONS FOR USE 

TO RUN THE SIMULATION WHICH VARIES THE NUMBER OF INSTRUCTIONS 
BETWEEN EACH NODE. 

A. NORMAL DATA INCLUDES: 


A. NORMAL DATA 


(1) . SIMNUM 

(2) . MINPJT 

(3) . MEANLN 

(4) . NUMOUT 

(5) . IMREAD 

ARC IS 
INREAD 

(6) . OELADD 

(7) . ITIME = 

(8) . MMTTR = 30 

(9) . MEANER = 20 

(10) . MEANIT = 20 

(11) . AVCHNG = 3 

( 12) . ASTER = ****•• 
(13) . N = 30 


INCLUDES: 
= 2 
= 1 
= 10 
= 30 
= 0 
TO BE 
- 1 
= 0 
50 


IF THE NUMBER OF INSTRUCTIONS IN EACH 
DETERMINED RANDOMLY - OTHERWISE 


IN MILLISECONDS) 
IN MINUTES) 


c 

(1). SIMNUM 

= 1 

C 

c 

( 2) . MINPUT 

= 3 

c 

c 

(3). MEANLN 

= 1 

c 

c 

(4). NUMOUT 

= 30 

c 

c 

(5). INREAD 

= 0 

c 

c 

(6). DELADD 

= 0 

c 

c 

17). ITIME = 

50 (IN MILLISECONDS) 

c 

c 

(8) . MMTTR * 

30 (IN MINUTES) 

c 

c 

(9). MEANER 

= 20 

c 

c 

(10). MEANIT 

= 20 

c 

c 

(11). AVCHNG 

= 3 

c 

c 

(12). ASTER = 


c 

c 

(13). N = 30 


r 

w 

c 

B. INPUT ADJACENCY MATRIX 

c 

c 



c 

c 

2. TO RUN THE SIMULATION WHICH VARIES THE NUMBER OF 

INPUTS USED: C 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OOOO 








noooononoooooooooonoooononooononoooooooonooooooo 


B. INPUT ADJACENCY MATRIX 

C. IF INREAD = 1 INPUT MATRIX OF ARC LENGTHS. 

3. TO RUN THE SIMULATION WHICH VARIES THE NUMBER OF ARCS BETWEEN 
THE NODES. 

A. NORMAL DATA INCLUDES: 

( 1) . SIMNUM = 3 

(2) . MINPUT = 3 

(3) . MEANLN = 10 

(4) . NUMGUT = 20 

(5) . INREAD = 0 

(6) . DELADD =0 IF DELETING ARCS FROM STRUCTURE. 

DELADD = 1 IF ADDING ARCS TO THE STRUCTURE. 

(7) . ITIME = 50 (IN MILLISECONDS) 

(8) . MMTTR = 30 (IN MINUTES) 

(9) . MEANER = 20 

(10) . MEANIT = 20 

(11) . AVCHNG = 3 

(12) . ASTER = ***** 

(13) . N = 30 

B. IF DELADD = 0 INPUT MOST COMPLEX STRUCTURED ADJACENCY 
MATRIX. 

IF DELADD = 1 INPUT SIMPLEST STRUCTURED ADJACENCY MATRIX. 

C. INPUT 20 (NUMOUT) PAIRS OF VALUES N1 AND N2 WHICH 
DETERMINE THE ARC TO BE DELETED OR ADDED. 

4. TO RUN THE SIMULATION WHICH VARIES THE NUMBER OF LOOPS: 

A. NORMAL DATA INCLUDES: 

(1) . SIMNUM = 4 

(2) . MINPUT = 3 

(3) . MEANLN = 10 

(4) . NUMOUT = 20 

(5) . INREAD = 0 

(6) . DELADD =0 IF DELETING LOOPS FROM STRUCTURE. 

DELADD =1 IF ADDING LOOPS TO THE STRUCTURE. 

(7) . ITIME = 50 (IN MILLISECONDS) 

(8) . MMTTR = 30 (IN MINUTES) 

( 9) . MEANER = 20 

(10) . MEANIT = 20 

(11) . AVCHNG = 3 

(12) . ASTER = ***** 

( 13). N = 30 

B. IF DELADD = 0 INPUT MOST COMPLEX STRUCTURED ADJACENCY 
MATRIX. 

IF DELADD = 1 INPUT SIMPLEST STRUCTURED ADJACENCY MATRIX. 

C. INPUT 20 (NUMOUT) PAIRS OF VALUES N1 AND N2 WHICH 
DETERMINE THE LOOP TO BE DELETED OR ADDED. 


oooo 




oooooooooooooooonoooooo 


* 


TO 

A. 


RUN THE SIMULATION WHICH VARIES THE SEED TO THE INPUT 
NORMAL DATA INCLUDES: 

*. 5 


( I) 
(2) 

(3) 

(4) 
( 5) 


B. 

C. 

D. 


(6) 

(7) 

(8) 
(9) 

( 10 ) 
( 11 ) 
( 12 ) 
(13) 
INPUT 


SIMNUM 

MINPUT 

MEANLN 

NUMOUT 

INREAD 

ARC IS 

INREAD 

DELADD 

I TI ME = 

MMTTR = 


1 

= 10 
= 30 
= 0 
TO BE 
= 1 
= 0 
50 
30 


MEANER = 20 
MEANIT = 20 
AVCHNG = 3 
ASTER = ***** 

N = 30 

ADJACENCY MATRIX 


IF THE NUMBER OF INSTRUCTIONS IN EACH 
DETERMINED RANDOMLY - OTHERWISE 

(IN MILLISECONDS) 

(IN MINUTES) 


IF INREAD = 1 
IZ = 1722632 


NOTE: LAST DATA CARD 


,^^Ja^v« X ust 4 ?e 6D le ^ T s HS on previous run. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


INUIC. LWJ I UK I H Q 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMMON/ ALLL/N# NODES (30 f 30 ) ACTn „ mm _ AwrwMG 


AVCHNG 

INST* 


100 

108 

140 

151 


CCMMON/ARRAY/X(30*30)* NUMPTS(30), ASTQUT(30)* 
COMMON/ERROR/ISEED(30*30)* ITER(30,30)* MEANER* 

CGMMC N/OUT/NSEED » INPUT, SVSEED(30,30) 
CCMMON/GEN/MEANLN* IX 
COMMON/NEW/IW 

CCMMON/GR AF/ INREAD TIIU „, sn , Tu-ntiMTl^n* 

DIMENSION TTIME(50), TREP(50)* TLINK(50)» TKOUNTl30) 

DIMENSION PSEED(50 )» FJ5§)t FNA(50)» FLA(50 ) 

DIMENSION TSEED( 50 ), TINST(_50), PTEST(50) 

INTEGER*4 CHANGE* 

INTEGER*4 DELADD 
REAL*4 MTTR 
(13 ) 

( 110 ) 

(1615) 

I'l't 13X * 

EACH NODE 


MEANIT, IZ 


TINPUT(50) 


AVCHNG, TESTED(30,30), TEST, SIMNUM, SVSEED 


( ' 1 


13X 


‘THIS 

•THIS 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
1E TWEEN 

152 FORMAT 
1S ED. 1 ) 

153 FORMAT 
1 HE NODES.*-) 

154 FORMAT (•1', 

155 i F determines’THE"INPUT'PATH.•) 


RUN WILL VARY THE NUMBER OF INSTRUCTIONS B 

RUN WILL VARY THE NUMBER OF INPUTS TO BE U 

{•1', 13X, ’THIS RUN WILL VARY THE NUMBER OF ARCS BETWEEN T 


13X 

13X, 


THIS 

•THIS 


RUN 

RUN 


WILL 

WILL 


VARY 

VARY 


THE 

THE 


NUMBER OF LOOPS.') 

RANDOM NUMBER SEED WHICH 


i 





156 FORMAT < * 0•, 
1 READ IN AND 
180 FORMAT ( 'O' , 
203 FORMAT ( 'O' , 
1 ZERO*//) 

“ ('O' 


13X, * THE NUMBER 
IS NOT RANDOM.•) 


OF INSTRUCTIONS IN EACH ARC HAS BEEN 


13X» ‘ADJACENCY MATRIX'//7X» 30(14)) 

13X» 'NOTE: NONEXISTENCE 3F AN ARC IS INDICATED BY A 


, 13X, 'THERE WILL BE', 13, * GRAPHS EVALUATED THIS RUN 


204 FORMAT 

1 .*//) 

205 FORMAT 

206 FORMAT 

207 FORMAT 
IN TIME 

2 OF', 13,' MILLISECONDS.'//13X, 'MEAN TIME 

3 AN EXPONENTIAL DISTRIBUTION WITH A MEAN OF', 13, 

4* THE NUMBER OF INPUT PATHS BEING INVESTIGATED IS', 13, 
5'ERRORS ARE SEEDED BY AN EXPONENTIAL DISTRIBUTION WITH 


(7(15, 5X)) 
(A5) 

('O', 12X, 
PER 


•INPUT ASSUMPTIONS (VARIABLES):'//13X, 'EXECUTIO 
INSTRUCTION HAS AN EXPONENTIAL DISTRIBUTION WITH A MEAN 


6 13,' INSTRUCTIONS PER 
7CTI0NS CHANGED WHEN AN 
208 FORMAT ('O', 13X, 'THE 
1D ISTRIBUTED FROM 1 TO* 


TO REPAIR AN ERROR HAS 
' MINUTES.'//13X, 
•.'//13X, 

A MEAN OF* 


ERROR.•//13X, 'THE AVERAGE NUMBER OF INSTRU 
ERROR IS FOUND IS', 13,'.*) 

NUMBER OF ITERATIONS PER LOOP IS UNIFORMLY 

13, '.*) 


209 FORMAT ('O', 13X, 'THE MEAN NUMBER OF INSTRUCTIONS BFTWEEN EACH NO 


IDE IS 
227 FORMA 
230 FORMAT( 


EXPONENTIALLY DISTRIBUTED WITH 
('O', 13X, 'SAMPLE INPUT PATH 
•O', 11X, 14, 4X, 14) 


A MEAN OF', 
FROM*//13X, 


13, '.' ) 

•NODE TO NODE'//) 


13X, 


240 FORMAT('O' 

1 F 8.2, 

2 // 13X, 'AND AN 
245 FORMAT ('O', 13X, 
250 FORMAT! 

lSEEDED 
255 FORMAT! 

260 FORMAT( 

265 FORMAT( 


FOR 

c 


NUMBER*,13,• 


THE EXECUTION TIME 
TIME OF*, F8.2, • 
TIME OF' , F8.2, • 


INPUT 

WITH A TOTAL REPAIR 
AVERAGE LINK TRAVERSAL 
•ERROR FOUND'//) 

I', 13X, 14, • SEEDED ERRORS REMAINING'/■/13X, 

ERRORS REMAINING'//7X» 30(14)) 

• », 7X, 30(A4)) 

'O', 15. ' *' 

• 1 * - 


IS* , 

HOURS' 
SECONDS') 


•MATRIX OF 


11X. 


, 

13, 


, 30(14)) 
INPUTS'//13X, 


13 




*//13X, F6.2, ' 


113, ' ERRORS SEEDED'//13X, 

2 • PERCENT RESIDUAL ERRORS 

3 //) 

O', 9X, F10.5, * SECONDS 
HOURS TO REPAIR ERRORS' 
TRAVERSAL TIME'//) 


14, • INSTRUCTIGNS'//13X, 


RESIDUAL 


ERRORS'//13X, F6.2, 
PERCENT OF ARCS TESTED* 


2 70 FORMAT ( 
1 F10.6, 
2RAGE ARC 
275 FORMAT( 

1 OF ARCS 

2 F5.2, ' 

280 FORMAT ( 

281 FORMAT ( 

282 FORMAT ( 

283 FORMAT ( 

284 FORMAT ( 

285 FORMAT ( 

286 FORMAT ( 


EXECUTION TIME'//12X, 

//12X, F10.6, • SECONDS AVE 


u-, 12X, F5.2, 
// 13X, F5.2, * 
IS THE RATIO OF 


0'//13X , 
0'//13X , 
0'//13X , 
0 *//13X » 
1«, IX) 
0'//l3X, 
0 *//13X , 


IS THE RATIO OF ACTUAL TO MAXIMUM NUMBER 
IS THE RATIO OF NODES TO ARCS'//13X, 

LOOPS TO ARCS'//) 


•HISTOGRAM FOR 
'HISTOGRAM FOP 
•HISTOGRAM FOR 
•HISTOGRAM FOR 

EXECUTION TIME 
EXECUTION TIME 


EXECUTION TIME') 

REPAIR TIME’) 

PERCENT RESIDUAL ERRORS' 
PERCENT ARCS TESTED') 


VS RESIDUAL ERRORS') 

VS PERCENT RESIDUAL ERRORS') 









oo 


I 


« 


ON 


287 FORMAT 

288 FORMAT 

289 FORMAT 

290 FORMAT 
1SIDUAL 

291 FORMAT 

292 FORMAT 

293 FORMAT 
1BER OF 

294 FORMAT 

295 FORMAT 

296 FORMAT 

297 FORMAT 
1RCENT 

298 FORMAT 
ID' ) 

299 FORMAT 

ID* ) 

300 FORMAT 

301 FORMAT 

302 FORMAT 

303 FORMAT 


< *0'//13X, 
<'0'//13X, 
('0'//13X * 
('0'//13X , 
ERRORS*) 

<'O'//13X t 
(*0'//13X, 
( '0V/13X, 
ARCS' ) 

('0 *//13X , 
(*0'//13X, 
CO ' // 13X , 
( ' 0'//13X , 


•REPAIR TIME VS RESIDUAL ERRORS') 

•NUMBER OF INPUTS VS RESIDUAL ERRORS') 

•NUMBER OF INSTRUCTIONS VS RESIDUAL ERRORS') 
•RATIO OF ACTUAL TO MAXIMUM NUMBER OF ARCS VS 


RE 


•RATIO OF 
'RATIO OF 
EXECUTION 


NODES TO ARCS VS 
LOOPS TO ARCS VS 
TIME VS RATIO OF 


RESIDUAL ERRORS *l 
RESIDUAL ERRORS*) 
ACTUAL TO MAXIMUM 


NUM 


VS 


OF ARCS TESTED' ) 

(* O'//13X , 'RATIO 


EXECUTION TIME 
EXECUTION TIME 
•NUMBER OF INPUTS 
•RATIO OF ACTUAL 


RATIO OF NODES TO ARCS') 

RATIO OF LOOPS TO ARCS') 

VS PERCENT OF ARCS TESTED*) 

TO MAXIMUM NUMBER OF ARCS VS PE 


OF NODES TO ARCS VS PERCENT OF ARCS TESTE 


('0'//13X» 'RATIO OF LOOPS TO ARCS VS PERCENT OF ARCS TESTE 


1 

2 

3 

4 

5 

6 
7 


( 1 0•//13X * 'NUMBER OF INPUTS VS 

C0'//13X, 'NUMBER OF INPUTS VS _ _ 

(•0'//13X»'EXECUTION TIME VS NUMBER OF INSTRUCTIONS') 


EXECUTION TIME') 

PERCENT RESIDUAL ERRORS') 


•DATA SUMMARY 

3X, 'REPAIR' 

< 


'//62X, 'PERCENT' 
4X» 'ACTUALS 4X, 
’ * ™- 'RESIDUAL' 


ERRORS', 3X, 

, 6X , 'TIMES 5X, 'TO MAX' 


IS 59X, 

2X» •EXECUTION' , 

' LOOPS'/4X* 'RUN*, 35X, 

4X t 'ARCSS 6X, 'TIME* 

•T0*/2Xt 

5X, 'NODES', 4X* 'INPUTS’, 4X , 

'- 4X, 'ERRORS', 4X, 'TESTED*, 4X, 
•ARCS', 6X, 'ARCS', 6X, 'ARCS'/) 

« *« -- - -- - 


3X 

'NODES’, 5X, 


'PERCENT' 

i 


2X» 
6X, 


RESIDUAL 1 


TO' 


8X i 


•NUMBER' 

« 


4X, 'ERRORS 1 , 


, 4X, 
(SEC)' 


INSTR. 

_ _ _ , .... _ , 4X, • 

( HRS)' , 4X, 

304 FORMAT ('O', 2X, 13, IX, •*', 5X, 13, 7X, F5.0, 3(5X, 

1 2(4X, F6.2), 2X, 2F10.5, 3X, F5.3, 2(5X, F5.3)) 

305 FORMAT ('O', 13X, 'THE LAST SEED USED WAS', Ill, '.') 
READ IN THE TYPE OF SIMULATION TO BE RUN (SIMNUM) 

(5,100) SIMNUM 

IN NUMBER OF INPUTS TO BE USED (MINPUT) 

(5,100) MINPUT 

THE MEAN NUMBER OF INSTRUCTIONS BETWEEN EACH 
,100) MEANLN 

THE NUMBER OF GRAPHS TO BE EVALUATED (NUMOUT) 
,100) NUMOUT 

WHETHER OR NOT THE NUMBER OF INSTRUCTIONS IS TO 
,100) INREAD 

WHETHER DELETING AN ARC FROM THE STRUCTURE 
ADDING AN ARC TO THE STRUCTURE. 

READ (5,100) DELADD 

READ IN THE MEAN TIME TO EXECUTE AN INSTRUCTION IN 
(ITIME) 

READ (5,100) ITIME 


SEEDED*, 
7X, 


F5.0), 


READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 


IN 

(5 

IN 

(5 

IN 

(5. 

IN 


NODE (MEANLN) 


BE READ IN 


OR 


MILLISECONDS 





no 


\ 




57 

59 


61 

62 


READ IN MEAN 
READ (5,100) 

READ IN MEAN 
READ (5,100) 

READ IN MEAN 
READ (5,100) 

READ IN MEAN 
FOUND (AVCHNG) 

READ (5,100) AVCHNG 
READ IN THE ASTERICKS 
READ (5,206) ASTER 
READ IN THE NUMBER OF NODES (N) 
READ (5,100) N 
READ IN THE ADJACENCY 
DC 20 1 = 1,N 

5,140) (NODES ( I 


TIME TO REPAIR AN ERROR IN MIN (MMTTR) 

MMTTR 

NUMBER OF INSTRUCTIONS THAT ARE ERROR FREE (MEANER) 
MEANER 

NUMBER OF ITERATIONS FOR EACH LOOP (MEANIT) 

MEANIT 

NUMBER OF INSTRUCTIONS CHANGED WHEN AN ERROR IS 


TO PUT AROUND OUTPUT ARRAYS (ASTER) 


MATRIX (NODES) 


20 READ ( 


IF (SIMNUM.EQ.il WRITE 

IF (SIMNUM.EQ.2) WRITE 

IF (SIMNUM.EQ.3) WRITE 

IF (SIMNUM.EQ.4) WRITE 

IF (SIMNUM.EQ.5) WRITE 

IF (INREAD.EQ.l) 

WRITE (6,207) ITIME, 
WRITE (6,208) MEANIT 
WRITE (6,209) MEANLN 
WRITE (6,204) NUMOUT 
NNOUT = 1 
IX = 71286223 
IW = IX 

DETERMINE THE MAXIMUM 
MAXARC = N*(N-1) 

NUMPTS IS THE LABEL FOR 
1, N 
= I 

ASTOUT ARE TO 
= ASTER 


J) 


(6 
(6 
(6 
(6 
( 6 


J=1 , 

151 ) 

152 ) 

153) 

154) 
155 ) 


N) 


WRITE(6,156) 

MMTTR, MINPUT, 


MEANER, AVCHNG 


DO 57 I 
NUMPT S(I) 

ASTER AND 
ASTOUTII ) 
CONTINUE 
CALL DIGRAF 
DO 62 I = 1, 

DC 61 J = 1, 
ISEED IS ~ “ 

I S E ED(I , J 
CONTINUE 
CONTINUE 
SEED THE 
IF (SIMNUM.EQ.l) 
CALL SEED 
THIS SEED CAN BE 


NUMBER OF ARCS IN THE ADJACENCY MATRIX 
THE OUTPUT ARRAY 


PUT ASTERICKS AROUND THE ARRAY 


N 
N 

THE ARRAY 
) = 0 


OF ERRORS SEEDED 


PROGRAM 


WITH 
IX = 


ERRORS 

1722632 


CHANGED TO THE LAST VALUE OF THE SEED ON THE 
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C 

C 

C 

c 

c 

c 

c 


c 


c 


63 


A CONTINUATION 
READ (5,108) 


OF 

IX 


THE DATA IS DESIRED 


I 

NUMBER 
X( I,J) 
NUMBER 


GO TO 65 
OF ARCS IN 


THE PROGRAM 


OF INSTRUCTIONS 


65 


68 


77 


78 


79 


80 


82 

83 


OF 

+■ ISEEDU 


ERRORS 

t J) 


SEEDED 


).EQ.l) LOOPS = LOOPS + 1 


PREVIOUS TEST IF 
IF ( SIMNUM.EQ.5) 

ARCS = 0. 

INST = 0 
NSEED = 0 
DC 65 1=1,N 
DO 65 J = 1,N 

ZERO THE MATRIX WHICH RECORDS WHETHER OR NOT AN ARC HAS BEEN TEST 
iESTED( I , J) = 0 
IF (NODESt I,J).NE.1) 

CALCULATE THE NUMBER 
ARCS = ARCS + 

CALCULATE THE 
INST = INST * 

CALCULATE THE 
NSEED = NSEED 
CONTINUE 

CALCULATE THE NUMBER OF LOOPS IN THE PROGRAM 
LCGPS = 0 
DC 68 I = 1,N 
DO 68 J = 1 •r 
IF ( NODES ( I, J 
CONTINUE 

CHOOSE SAMPLE INPUT 
INPUT = 1 
TCTREP = 0. 

TIME = 0. 

AVLINK = 0. 

TTI ME(NNQUT) = 0. 

TREP(NNOUT) = 0. 

TLINK(NNOUT) = 0. 

IZ = IX 
NCDE = 1 
WRITE (6,227) 

NUMSUC IS THE NUMBER OF 
NUMSUC = 0 
DO 80 J = 1 » N 

IF (NGDESiNODE,J).EQ.l) 

CONTINUE 

IF (NUMSUC.EQ.O) GO TO 
IF (NUMSUC.NE.l) GO TO 
K = 1 
GO TO 83 

THE SUCCESSORS ARE CHOSEN RANDOMLY WITH A UNIFORM DISTRIBUTION 
CALL RANDOMtIZ, U, 1) 

K = 1 + NUM5UC*U 
KK = 0 

DO 85 J = 1, N 


SUCCESSORS 


NUMSUC = NUMSUC + 1 


96 

82 






on o o o o oo 


r 


4 


KK = KK + 1 


XTIME/ 1.0E03 


PATH 
TO 95 


BE INSERTED 


IF (NODES!NODE,J}.EQ.l) 

IF (KK.EQ.K) GO TO 86 

85 CONTINUE 

86 L = J 

TIME IS ASSUMED TO HAVE A MEAN OF 
INSTRUCTION WITH AN EXPONENTIAL 
CALL EX PON( I Z» XTIME, 1) 

XTIME = XTIME * ITI ME 
IF (NODE.L E.L) GO TO 88 
CONVERT TIME INTO HOURS 
TIME = TIME * XTIME * ITER(NODE »L) 

1 1.0E03 
GO TO 89 

88 TIME = TIME + X(NODE,L) * 

89 WRITE (6,230) NUDE,L 
TESTED!NODE,L) = I 
IS THERE AN ERROR IN THIS 
IF!ISEED(NODE,L).EQ.O) GO 
ERROR FOUND 
WRITE (6,245) 

SEE IF A NEW ERROR SHOULD 
CALL NU SEED 

90 I SEED!NODE,L) = I SEED(NODE,L) - 
CALL E X PON(IZ, U, 1) 

MEAN TIME TO REPAIR HAS AN 

MEAN OF MMTTR MINUTES 
MTTR = MMTTR * U 
TOTREP IS THE TOTAL REPAIR 
TCTREP = TOTRCP + MTTR/ 60. 

GO TO 78 

95 NODE = L 
GO TO 79 

AVLINK IS THE AVERAGE ARC 

96 AVLINK = (TIME/ARCS) 

WRITE (6,240) INPUT, TIME, 

IF (SIMNUM.EQ.2) GO TO 97 
TTI ME(NNGUT) = TTIME!NNOUT) 
TREP(NNOUT) = TREP(NNOUT) + 
TLINK(NNOUT) = T LINK!NNOUT) 

GO TO 99 

97 NOUT = NNOUT - 1 

IF (NNOUT.GT.l) GO TO 98 
TLINK(I) = AVLINK 
TREP(l) = TOTREP 
TTIME(l) = TIME 
GO TO 99 

98 TLINK(NNOUT) = TLINK(NOUT) + AVLINK 
TREP(NNOUT) = TREP(NOUT) + TOTREP 


ITIME MILLISECONDS PER 
DISTRIBUTION 


* (X (NODE, L) X( L,NODE ) ) / 


EXPONENTIAL DISTRIBUTION WITH A 


TIME IN HOURS 


TRAVERSAL TIME IN MINUTES 


TOTREP, AVLINK 

+ TIME 
TOTREP 
+ AVLINK 
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t 


4 


I 


TTIME(NNGUT) = TTIME(NOUT) + TIME 
99 TIME = 0. 

TCTREP = 0. 

AVLINK = 0. 

INPUT = INPUT + 1 

C KCUNT IS A COUNTER OF THE NUMBER OF SEEDED ERRORS REMAINING 

KCUNT = 0 
DO 101 I = 1 , N 
DO 101 J = 1 , N 

IF (NODES(I,J).EQ.l) KOUNT = KOUNT + ISEED(I,J) 

101 CONTINUE 

WRITE (6,250) KOUNT, (NUMPTS(I), I = 1,N) 

WRITE (6,255) (ASTOUT(I), I = 1,N) 

DO 110 I = 1t N 

110 WRITE (6,260) NUMPTS(I) , (ISEED(I,J), J=1,N) 

C DISTINCT SEED FOR RANDOM NUMBER GENERATOR FOR EACH INPUT 

IX = IX - 12345 
IF (SIMNUM.EQ.2) GO TO 111 
IF (SIMNUM.EQ.5) GO TO 111 

C MINPUT IS THE MAXIMUM NUM3ER OF INPUT PATHS TO BE CHECKED 

IF (INPUT.LE.MINPUT) GO TO 78 
C CALCULATE THE NUMBER OF ARCS TESTED 

111 TEST = 0 

DC 112 I = 1,N 
DO 112 J = 1,N 

IF (TESTED(I,J).EQ.l) TEST = TEST + 1 

112 CONTINUE 

INPUT = INPUT - 1 

C PTEST IS THE PERCENT OF ARCS TESTED 

PTEST(NNOUT) = 100. * TEST / ARCS 
C F IS THE RATIO OF ACTUAL TO MAXIMUM NUMBER OF ARCS 

F(NNOUT) = ARCS / MAXARC 
C FNA IS THE RATIO CF NODES TO ARCS 

FNA(NNOUT) = N / ARCS 
C FLA IS THE RATIO OF LOOPS TO ARCS 

FLA(NNOUT) = LOOPS / ARCS 
TINST(NNOUT) = INST 
TINPUT(NNOUT) = INPUT 
TSEEO(NNQUT) = NSEED 
TKOUNT(NNOUT) = KOUNT 

PSEED(NNOUT) = 100. * (TKOUNT(NNOUT)' / NSEED) 

WRITE (6,265) INPUT, INST, NSEED, KOUNT, PSEED(NNOUT), 

1 PTEST(NNOUT) 

WRITE (6,270) TTI ME(NNOUT), TRE P(NNOUT), TLINK(NNOUT) 

WRITE ( 6,275) F(NNOUT), FNA(NNOUT 5, FLA(NNOUT) 

NNOUT = NNOUT + 1 

IF (NNOUT.GT.NUMOUT) GO TO 118 

IF (SIMNUM.EQ.l) MEANLN = MEANLN + 1 
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* 


( 


IF (SIMNUM.EQ.l) GO TO 57 
IF (SIMNUM.EQ.2) INPUT = INPUT + I 
IF (SIMNUM.EQ.2) GO TO 78 
00 113 I = 1,N 
DC 113 J = 1,N 

113 ISEED(I,J ) = SVSEED(I,J) 

IF (SIMNUM.EQ.5) GO TO 63 

115 IF (DELADD.EQ.0) CALL DELARC 
IF (DELADD.EQ.1) CALL ADDARC 

116 WRITE (6,180) (NUMPTS(I), 1=1,N) 

WRITt (6,2551 (ASTOUT(I), I = 1,N) 

DO 117 I = 1 , N 

117 WRITE (6,260) NUMPTS(I), (NODES(I,J), J=1,N) 
WRITE (6,203) 

GO TO 63 

118 WRITE (6,305) IL 
WRITE (6,284) 

CALL PLOTP (TTIME, 

WRITE (6,285) 

WRITE (6,284) 

CALL PLOTP (TTIME, 

WRITE (6,236) 

WRITE (6,284) 

CALL PLOTP (TREP , 

WRITE (6,287) 

IF (SIMNUM.NE.2) 

WRITE (6,284) 

CALL PLOTP (TINPUT 
WRITE (6,288) 

WRITE (6,284) 

CALL PLOTP (TINPUT, 

WRITE (6,296) 

WRITE (6,284) 

CALL PLOTP (TINPUT 
WRITE (6,300) 

WRITE (6,284) 

CALL PLOTP (TINPUT 
WRITE (6,301) 

119 IF ((SIMNUM.EQ.2). 

WRITE (6,284) 

CALL PLOTP (TTIME, 

WRITE (6,302) 

WRITE (6,284) 

CALL PLOTP (TINST, 

WRITE (6,289) 

IF (SIMNUM.EQ.l) 

WRITE (6,284) 

CALL PLOTP (F, 


TKOUNT, 

NUMOUT, 0) 

PSEED, 

NUMOUT, 0) 

TKOUNT, 

NUMOUT, 0) 

GO TO 119 

, TKOUNT, 

NUMOUT, 0)- 

, PTEST, 

NUMOUT, 0) 

, TTIME, 

NUMOUT, 0) 

, PSEED, 

NUMOUT, 0) 

OR.(SIMNUM 

I.EQ.5)) GO 

TINST, 

NUMOUT, 0) 

TKOUNT, 

NUMOUT, 0) 

GO TO 121 

TKOUNT, 

NUMOUT, 0) 








122 


WRITE (6,290) 

WRITE (6,284) 

CALL PL OTP (FNA, 
WRITE (6,291) 

WRITE (6,284) 

CALL PLOTP (FLA, 
WRITE (6,292) 

WRITE (6,284) 

CALL PLOTP (TTIME, 
WRITE (6,293) 

WRITE (6,284) 

CALL PLOTP (TTIME, 
WRITE (6,294) 

WRITE (6,234) 

CALL PLOTP (TTIME, 
WRITE (6,295) 

WRITE (6,234) 


TKOUNT, 

NUMOUT, 

0) 

TKOUNT, 

NUMOUT, 

0) 

F, 

NUMOUT, 

0) 

FNA, 

NUMOUT, 

.0) 

FLA, 

NUMOUT, 

0) 


CALL PLOTP (F, 

WRITE (6,297) 

WRITE (6,284) 

PTEST, 

NUMOUT, 

0) 

CALL PLOTP (FNA, 

WRITE (6,298) 

WRITE (6,284) 

PTEST, 

NUMOUT, 

0) 

CALL PLOTP (FLA, 

WRITE (6,299) 

DO 125 K = 1,2 

WRITE (6,303) 

PTEST, 

NUMOUT, 

0) 

WRITE (6,255) (ASTOJT(I), I 
DO 122 1 = 1,NUMOUT 

= 1 , N) 



121 


122 WRITE (6,304) I, N, TINPUT(I), TINST(I), TSEED(I), TKOUNT(I), 
1 PSFED(I), PTEST(I), TTIME(I), TREP(I), F(I), FNA(I), FLA(I) 
125 CONTINUE 

IF ((SIMNUM.NE.3). AND.. (SIMNUM.NE.4)) CALL GRAPHO 

NEED A MINIMUM OF 10 GRAPHS TO ANALYZE BEFORE CALLING HISTG 

IF (SIMNUM.NE.5) GO TO 126 

CALL HISTG (TTIME, NUMOUT, 0) 

WRITE (6,280) 

0 ) 


CALL HISTG (TREP, 
WRITE (6,281) 

CALL HISTG (PSEED, 
WRITE (6,282) 

CALL HISTG (PTEST, 
WRITE (6,283) 

126 CONTINUE 

STOP 


NUMOUT, 


NUMOUT, 0) 
NUMOUT, 0) 


END 





^ ^ I 
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ro 


SUBROUTINE DIGRAF 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C CONSTRUCT THE DIRECTED GRAPH 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COMMON/ALLL/N, N0DES(30,30) 

COMMON/ARR AY/X(■ 30 » 30) t NUMPTS(30)» AST0UT(30), AVCHNG 
CCMMON/GEN/MEANLN, IX 

C0MM0N/ERR0R/ISEED(30»30), ITER(30,30), MEANER, INST, MEANIT, IZ 
C OMMON/GRAF/INREAD 


100 
130 
140 
150 
160 
1 80 
190 
200 


INTEGERS CHANGE, AVCHNG 
FORMAT (16F5.0) 


14, IX, 'NODES IN DIRECTED GRAPH.'//) 
•MATRIX OF ARC LENGTHS•//7X, 30(14)) 

15, ' *' , 30(F4.0)) 

'NOTE: 0.0 INDICATES NO ARC FROM I TO 
'ADJACENCY MATRIX'//7X, 30(14)) 

15, • *», 30(14)) 

•NOTE: NONEXISTENCE 


C 


222 

225 

230 


10 


20 


O' 

O' 


7X, 
7X, 

13X 


30(A4)) 

30(A4)) 

15, • *• , 30(14)) 

•LOOP ITERATI0NS'//7X, 30(14)) 
RANDOMLY GENERATE INSTRUCTION LENGTHS 
READ IN INSTRUCTION LENGTHS 
GO TO 20 


FORMAT ('O', 10X. 

FORMAT('0' , 13X, 

FORMAT ('O', 

FORMAT ('O', 13X 
FORMAT ('O', 13X 
FORMAT ('O', 

FORMAT ('O', 13X 
1 ZERO'//) 

210 FORMAT (• 

FORMAT ( 

FORMAT (• 

FORMAT (•_ . . 

IF INREAD = 0 
IF INREAD = 1 
IF (INREAD.EQ.Oi 
DO 10 I = 1,N 
READ (5,100) (X(I 
GO TO 25 

GENERATE THE NUMBER OF INSTRUCTIONS 
EXPONENTIALLY DISTRIBUTED 
DO 25 1=1,N 
DO 25 J = 1,N 
X ( I , J ) = 0.0 

IF (NOOESt I,J).EQ.O) GO TO 25 
CALL EXPON (IX,XLNTH, 1) 

MEANLN IS THE EXPONENTIALLY DISTRIBUTED 
X ( I , J) = MEANLN * XLNTH + 1.0 
25 CONTINUE 

WRITE (6,130) N 

WRITE (6,140) (NUMPTS(I), I = 1,N) 

(AST OUT(I), I = 1,N) 


J'//) 


OF AN ARC IS INDICATED BY A 


J) , J = 1,N) 


BETWEEN EACH NODE - 


MEAN LENGTH OF EACH ARC 


WRITE { 

DO 30 
WRITE 
WRITE 
WRITE 


6,210 ) 
1=1 ,N 
(6,150) 
( 6 , 160 ) 
(6,180) 


NUMPT S(I), (X(I ,J) , J = 1,N) 
(NUMPTS(I), I =1 , N ) 


30 


ooo 







noon ooooo o oo 






* 


* 


WRITE (6,210) (ASTOUT(I), I = 1,N) 

DO 50 1=1,N 

50 WRITE (6,190) NUMPTS(I), <N0DES(I,J), J=1»N) 

WRITE (6,200) 

THE NUMBER OF ITERATIONS PER LOOP IS DISTRIBUTED WITH A UNIFORM 
DISTRIBUTION WITH A MEAN OF MEANIT 
DO 55 I = 1,N 
DO 55 J = 1,N 
I TER(I,J) = 0 

IF (NODESCI,J).EQ.0) GO TO 55 
I TER(I,J ) = 1 
IF (J.GT.I) GO TO 55 
CALL RANDOM (IX, TITER, 1) 

ITER(I,J) = 1 + MEAN IT*TITER 
55 CONTINUE 

CREATE LABELS FOR THE OUTPUT ARRAY 
WRITE (6,230) (NUMPTS(I), 1=1,N) 

WRITE (6,222) (ASTOUT(I), 1=1,N) 

DO 73 I = 1,N 

73 WRITE ( 6, 225) NUMPTS(I), (ITER(I , J), J=1,N) 

RETURN 

END 


SUBROUTINE ADDARC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

ADDARC ADDS AN ARC BETWEEN TWO NODES 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
COMMON/ALLL/N, NODES(30,30) 

COMMON/ARRAY/X(30,30), NUMPTS(30), ASTOUT(30)» AVCHNG 
COMMON/GEN/MEANLN, IX 
100 FORMAT (215) 

110 FORMAT I'l't 13X, 'A NEW ARC HAS BEEN ADDED FROM NODE', 13, 

1 • TO NODE', 13, ' OF LENGTH', F4.0) 

C ADD A NEW ARC FROM N1 TO N2 

READ (5,100) Nl, N2 
NODES(N1, N 2 ) = 1 
CALL EXPON (IX, XLNTH, 1J 

C MEANLN IS THE MEAN LENGTH OF THE ARC TO BE EXPONENTIALLY DISTRIB. 

X(M1,N2) = MEANLN * XLNTH +1.0 
WRITE (6,110) Nl, N2, X(N1,N2) 

RETURN 

END 


ooo 






t 


ro 

Vj\ 


SUBROUTINE DELARC 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 


c 

c 


DELARC DELETES AN ARC BETWEEN TWO NODES 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMMON/ALLL/N, N0DES(30,30) 

COMMON/ARRAY/XI30,30), NUMPTS(30), ASTOUT(30), AVCHNG 
C OMMON/ERROR/IS E E D(30 »3 0), ITER(30»30J, MEANER, INST, MEANIT, IZ 


I 


COMMON/OUT/NSEED 
INTEGER SVSEED 
100 FORMAT (215) 

110 FORMAT Cl', 13X 
• TO MODE*, 13, 

READ (5,100) Nl, 
NODES(N1.N2) = 0 
X(N1, N2J = 0.0 
ISEED(N1,N2) = 0 
SVSEED(N1, N2) = 
WRITE (6,110) Ml 
RETURN 
END 


INPUT, SVSEEDt30,30) 


N2 


AN ARC 
'/) 


WAS DELETED FROM NODE', 13, 


N2 


C 
C 

SUBROUTINE SEED 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C SEED THE GRAPH WITH ERRORS 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMMON/ALLL/N, N0DES(30,30) 

COMMON/ARRAY/X(30,30), NUMPTS(30), AST0UT(30), AVCHNG 

COMMON/ERROR/ISEED(30»30), ITER(30,30), MEANER, INST, MEANIT, IZ 

CGMMON/OUT/NSEED, INPUT, SVSEED(30,30) 

. C OMMON/GEN/MEANLN, IX 

INTEGER*4 CHANGE, AVCHNG, SVSEED 

IER(30) 

'SEED ERRORS AT INSTRUCTION INTERVALS'// 10X, 


DIMENSION ER(30) 

205 FORMAT (•O', 13X 
1 3014) 

210 FORMAT (•l•// 13X, 

220 FORMAT CO'// 13X, 

1'THE ERROR MATRIX«//7X, 


C 

C 


222 

225 


FORMAT 

FORMAT 


( 

(• 


O' 


7 X < 


'TOTAL NUMBER 
•TOTAL NUMBER 
30(14)) 


OF 

OF 


INSTRUCTIONS IS', 15) 

ERRORS SEEDED IS',14// 13X, 


30(A4)) 
15, • *• 


30(14)) 


GENERATE ERRORS WITH EXPONENTIAL DISTRIBUTION 
CALL EX PON(IX , ER, N) 


lOOOO OOOO 






non o non on 


SPREAD ERRORS TO GIVE MEAN OF MEANER 


DO 60 I = I»N 
IER(I) = MEANER * ER(I) 

60 CONTINUE 

NSEED 1 IS THE NUMBER OF ERRORS SEEDED 
NSEED = 0 

INST IS THE NUMBER OF INSTRUCTIONS 

NUMER - IS THE ERROR SEED/I NSTRUCTICN COMPARATOR 
NUMER = 0 
DO 70 I = 1,N 
DO 65 J = I * N 

IF (N0DES(I,J).EQ.O) GO TO 65 
COUNT THE NUMBER OF INSTRUCTIONS 

$Be e EB 8S M2 E T 5S°E§55i8 M & I? e JX?ions 

IF R ?J?LeIi) IERUI) = IER(II) / ( (XII,J) * ITER( I, J) )/MEANER) 

63 NUMER = NUMSP. + IER(II) 

IF (INST.LT.NUMER) GO TO 64 

I S EED( I ,J) = ISE ED(I,J) + 1 
NSEED = NSEED + 1 

II = II + 1 

IF (II.LE.N) GO TO 63 
I I = I 
GO TO 63 

64 NUMER = NUMER - IER(II> 

65 CONTINUE 

70 CONTINUE 

DO 71 I = 1 , N 

SVSEED J IS THE SAVE COPY OF THE ORIGINAL SEEDED ERRORS 

71 SV SE ED( I,J) = ISEEDtI,J) 

NSEED IS THE NUMBER OF ERRORS SEEDED 
WRITE (6,210) INST 

WRITE (6,205) (IER(I), I = 1,N) , 

WRITE (6,220) NSEED, (NUMPTS(I), I = 1,N) 

WRITE (6,222 ) f A STOUT(I)* 1=1,N) 


- IER(II> 


WRITE (6,220) NSEED, (NUMPTS(I), I = I,N) 
WRITE (6,222 ) f A STOUT(I)* 1=1,N) 

WRITE (6^225) NUMPTS(I), (ISEED(I,J), J=1,N) 


75 CONTINUE 


RETURN 
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SUBROUTINE NUSEED 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C CORRECTING PREVIOUS ERROR MAY CREATE A NEW ERROR C 

C C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCXCCCCCCCCCC 
COMMON/ALLL/N, NODES(30,30) 

COMMON/ARRAY/XI30,30), NUMPTS(30), AST0UTC30), AVCHNG 

COMMON/ERROR/ISEED(30,30), ITER(30,30), MEANER, INST, MEANIT, IZ 

COMMON/OUT/NSEED, INPUT, SVSEED!30,30) 

COMMON/NEW/IW 
INTEGER*4 CHANGE, AVCHNG 
DIMENSION ERR(2) 

247 FORMAT (*0‘, I3X, 'NEW ERROR INSERTED FROM NODE', 14,» TO NODE*, 

1 14//) 

NUM = 0 

C THE ERROR IS EQUALLY LIKELY TO RESULT IN CHANGES TO ALL ARCS 

CALL RANDOM (IW, ERR, 2) 

1=1+ N*ERR<1) 

J = 1 + N*ERR(2) 

C NEW ERROR CREATED AT (I,J) 

C IS THIS A NODE OF THE OIRECTED GRAPH 

10 IF (NODES!I,J).EQ.1) GO TO 20 
J = J + 1 
IF (J.GT.N) J = 1 
NUM = NUM + 1 

C IF NODE I IS A TERMINAL NODE NO NEW ERROR IS INSERTED 

IF (NUM.GT.N) RETURN 
GO TO 10 

C CHANCE OF INSERTING A NEW ERROR IS BASED ON THE NUMBER OF 

C INSTRUCTIONS THAT MUST BE CHANGED 

20 CALL RANDOM (IW, ERNU, 1) 

C AVCHNG IS THE ASSUMED (INPUT) AVERAGE NUMBER OF INSTRUCTIONS 

C CHANGED WHEN AN ERROR IS FOUND 

C CHANGE IS THE UNIFORMLY DISTRIBUTED NUMBER OF INSTRUCTIONS CHANGED 

. CHANGE = 1 + AVCHNG*ERNU*2, 

C CRATIO IS THE RATIO OF INSTRUCTIONS CHANGED TO THE NUMBER OF 

C INSTRUCTIONS IN THE UNIFORMLY RANDOMLY CHOSEN ARC 

CRATIO = CHANGE / X(I,J) 

C THE CRITERIA FOR INSERTING A NEW ERROR IS UNIFORMLY DISTRIBUTED 

CALL RANDOM (IW, RATIO, 1) 

C IF THE CHANGE RATIO (CRATIO) IS GREATER THAN THE CRITERIA (RATIO), 

C THEN A NEW ERROR IS INSERTED 

IF (CRATIO.LE.RATIO) RETURN 
50 I SEED( I , J ) = ISE ED ( I » J ) + 1 
NSEED = NSEED + 1 
WRITE (6,247) I,J 
RETURN 
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END 

SUBROUTINE GRAPHO 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C GRAPHICALLY DISPLAY THE DIRECTED GRAPH USING THE CALCOMP PLOTTER 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMMON/ALLL/Nt N0DES(30,30) 

DIMENSION ND'(30,30), X( 30), XX(2), Y(30), YY(2) 

DIMENSION XL(36) » YL(36), THETA(36) 

INTEGERS ITB U 2 J/ 12*0/ 

REAL-4 RTBI28)/28*0.0/ 

C ITB IS DRAWP ARRAY 

I T B U ) = 1 
I T 6 ( 2 ) = 0 
ITB(3) = 9 
IT8(4) = 15 
IT318) = 2 
ITB(91 = 0 
I TB{1C) = 2 
ITB{11) = 0 

C RTB IS DRAWP ARRAY 

RTB(l) = 1.0 
RTBI21 = 1.0 
DO 4 I = 1 y N 
DO 4 J = 1,N 

C ND IS ARRAY USED TO CHECK FOR PREDECESSORS 

ND(I» J) = 0 
4 CONTINUE 
XU) = 4.5 
Y(l) = 15.0 
^I 2 j — X(l) 

C CONST IS Y-DIRECTION SCALE FACTOR 

CONST = 30./N 
YI2) = Y(1) - CONST 
CALL DRAWP(2tXfY T ITB,RTB) 

XCGNST = 1.5 

C KOUNT IS X-DIRECTION SHRINK FACTOR - SHRINKS AFTER EVERY OTHER 

C FANOUT 

KCUNT = 1 

C LEVEL IS THE LEVEL IN THE DIRECTED GRAPH 

LEVEL = 2 
NN = N - 1 
ITB(1) = 2 

C LABEL IS THE NUMBER OF EACH NODE 

LABEL = 2 
DO 5 M = 1,N 

IF (NODES(1»M).EQ.O) GO TO 5 


noon 






1 


ro 

vo 


KK IS THE LAST NODE IN THE GIVEN LEVEL 
KK = M 
KT = KK + 1 

IF (NODES!l,KT).EQ.O) GO TO 6 

5 CONTINUE 

KM IS THE LAST NODE IN THE NEXT LEVEL 

6 DO 7 M = KK,N 

IF (NODES(KK,M).EQ.O) GO TO 7 
KM - M 
KX = KM + 1 

IF (NODES(KK»KX).EQ.O) GO TO 8 

7 CONTINUE 

8 DO 50 I = 2,NN 
11=1-1 

NUMSUC IS THE NUMBER OF SUCCESSORS 
NUMSUC =0 
DO 20 J = I,KM 

IF (NODES(I,J).EQ.OJ GO TO 20 
DO 18 K = 1,11 

IS THERE A NODE ALREADY IN EXISTENCE 
IF <ND(K,J).EQ.O) GO TO 18 

CONNECT NEW NODE TERMINATING AT EXISTING NODE 

17 XX(1i — X(I) 

XX<?) = X(J) 

YY(1) = Y(I) 

YY( 2) = Y( J) 

CALL DRAWP(2,XX,YY »ITB » RTB) 

ND(I,J) = 1 
GO TO 20 

18 CONTINUE 

NUMSUC = NUMSUC *■ 1 
20 CONTINUE 

IF (NUMSUC.LE.O) GO TO 44 
25 IF (MUMSUC.NE.U GO TO 30 
FANOUT TO ONE NEW NODE 
L = LABEL + 1 
LABEL = LABEL + 1 
ND(I,L) = 1 
X(L) = X( I ) 

Y(L) = Y(I) - CONST 
XX(1) = X(I) 

XX(2) = X(L) 

YY(1) = Y(I) 

YY(2) = Y(L) 

CALL DRAWP(2»XX,YY»ITB»RTB) 

IS THERE A NODE ABOVE THAT NEEDS CONNECTING TO THIS NODE 

DO 28 K = 1,11 

IF (NU(K,LI.EQ.l) GO TO 28 




IF (NODES(K»L).EQ.O) GO TO 28 
XX(I) = X(K) 

Y Y (1) = Y(K) 

XX( 2) = X(L) 

Y Y ( 2) = Y ( L ) 

ND(K,L) = 1 

CALL DRAWP(2»XX,YYtITBrRTB) 

28 CONTINUE 
GO TO 44 

FANOUT TO MULTIPLE NODES 

30 FAN = 2.*XCGN$T / (NUMSUC - 1) 

DO 35 L= L t NUMSUC 

M =' LABEL + L 
ND(I,M) = 1 
IF (L.GT.l ) GO TO 31 
X(M) = X(I) - XCONST 
GO TO 33 

31 LN = M - 1 

X(M) = X(LM) + FAN 

33 Y (M) = Y(I) - CONST 
XX(1) = X(I) 

YY< 1) = Y ( I ) 

XX{2) = X(M) 

Y Y { 2 ) = Y { M) 

CALL DRAWPI2»XX,YY»ITB,RTB) 

IS THERE A NODE ABOVE THAT NEEDS CONNECTING TO THIS NODE 
DO 34 K = 1* 11 

IF (ND(K,M).EQ.l) GO TO 34 
IF(NODES(K,M).EQ.Oi GO TO 34 
XXI1J = X(K) 

Y Y {1 ) = Y C K) 

XX(2) = X { M) 

YY( 2) = Y < M) 

ND(K*M) -- 1 

CALL DRAWP(2tXX,YY,ITB»RTB) 

34 CONTINUE 

35 CONTINUE 

LABEL = LABEL + NUMSUC 
K CUNT = KOUNT + 1 
40 IF (KOUNT.LT.2) GO TO 44 

REDUCE THE X-DIRECTION SCALE FACTOR 
XCCNST = XCONST/1.5 
KOUNT = 0 

IS THIS NODE STILL ON THE GIVEN LEVEL 
44 IF (I.LT.KK) GO TO 50 
NO* CHANGE LEVELS 
LEVEL = LEVEL + 1 

DETERMINE THE LAST NODE IN THE NEW LEVEL 
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c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


■EQ.O) GO TO 50 


DO 45 M = I,N 

IF (NODES! I,M) .EQ.O) GO TO 45 
KK = M 
KT = KK + 1 

45 1 ^0 F N!?^^uT“’ KT, • EQ • 0, G ° T ° 46 

, , DETERMINE the LAST NODE IN THE NEXT LEVEL AFTER THE PRESENT 

46 DO 47 M = KKjN 

IF (NODES(KK,M).EQ.O) GO TO 47 
KM = M 
KX = KM + I 
IF {NGDES,{KK?KX) 

47 CONTINUE 
50 CONTINUE 

GRAPH LOOPS 
DC 90 I = 2,N 

K = I - 1 

DO 80 J = 1,K 

IF (NODES! I,J).EQ.O) GO TO 80 

PI = 3.1415926 

ESTABLISH THE DIRECTION OF CURVE FOR LOOP 
IF (X(I).L T.X(J) ) PI = -PI 
STARTING POINT FOR CURVE 
XL(1) * XII! 

V L(1) = Y(I) 

ENDING POINT FOR THE CURVE 
XL(36) = X(J) 

Y L(3 6) = Y(J) 

vRic £ N 9,Y? IF 6? E ., THE DIFFERENCE between the X and Y VALUES 

AUlh - Al1J “ X ( J } 

YD IF = Y( I) - Y(J) 

R IS THE RADIUS OF THF LOOP 
R = 0.5 * SQRT!XDIF**2 f YDIF**2) 

(XC,YC) IS THE CENTER OF THE LOOP 
XC = (X( l ) + X(J))/2. 

YC = (Y( I ) + Y(J) )/2. 

THETE(l) IS THE INITIAL ANGLE 
THETA(l) = ATAN2(YDIF,XDIF) 

CREATE 36 POINTS IN A SEMICIRCLE FORMING THE LOOP 
DO 70 L - 2? 35 
M = L - 1 

THETA(L) = THETA(M) + PI/36. 

X L i L) = R * COS (THETA(L)J + 

TO * S ‘ N 1THETA,L »> 4 

80 ^^T L lS5E HP(36 ’ XL ’ ¥L ' 1TB ’ R ™ 

90 CONTINUE 


XC 

YC 


LEVEL 


RESP. 





C DRAW DIAMONDS AT EACH NODE 

IT8(I) = 3 

CALL DRAWP(N,X t Y,ITB,RTB) 

RETURN 

END 

//LINK.SYSPRINT DD SYS0UT=0 
//GO .FTC6FOOI DD $YS0UT = 0,DC8=(RECFM 
// UNIT=(SYSOUT,SEP=FT05FOOi),SPACE 
//GO.SYSIN DD * 

5 

1 

10 

50 

0 

1 

50 

30 

20 

20 

3 


30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 


1 0 

0 0 

0 1 

0 0 

0 0 

0 0 

0 0 

0 0 

0 1 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

1 o 

0 0 


0 0 
0 0 
1 0 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 

0 0 

0 0 

0 0 

1 o 

0 0 

0 1 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


FBA»LRECL=133» BLKSIZE=3325)t 
(CYLf (4« in 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
1 1 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 0 0 
0 0 0 
0 0 0 0 
0 0 0 
0 0 0 0 
0 0 0 
0 0 0 0 
0 0 0 
0 0 0 0 
0 0 0 
0 0 0 0 
0 0 0 
1 0 0 0 
0 0 0 


Olio 

0 0 0 



0 0 0 0 
0 0 0 
0 0 0 0 
0 0 0 
0 0 0 0 


0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 






ooooooooooooooooo 


OOOOOOOOOOOOOOOOO 


OOO OOOOOOOOOOOOOOOOO OOOOOOOOOOOO-^OO 


OOO OOOOOOOOOOOOOOOO00000—100000^0000 


OOOOOOOOOOOOOOOOOOhOOOOOOOOOhOOOOOO 


OOO OOOOOOO—'0000000000000-^0000000000 


OOO OOO O OOO OOOOO OOOOO 00-^0 OOO OOO OOOOO 


OOOOO OO0000000000000—1000 OOO OOO O OOO o 


o-^oooooooooooooo-<oooooooooooooooooo 


OOOOOOOOO—tooHOHOOOOOOOOOOOOOOOOOOOO 


OOOOOOOO—'OOOOOOOOOOOOOOOOOOOOOOOOOO 


OOO—*»—«0000000OOO OOOO OOOOOOOOOOOOOOOO 


OO—*000000000000OOOOOOOOOOOOOOOOOOOO 


—•OOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


133 


622150796 







LIST OF REFERENCES 




« 




# 


k 




1. Boehm, B. W. f " Software and Its Impact: A Quantitative 
Assessment' 1 , Dat amation, v. , p. 48-59, Bay 1973. 

2. Boehm, B. W., "The High Cost of Software", Proceedings 
of a S ympo si um on the High Cost of So ftwa re, v. Jack 
Goldberg (ed), Stanford Research Institute, p. 27-40, 
1973. 

3. Bradley, G. H., Green, T. F., Howard, G. T., and 

Schneidewind, N. f., "Structure and Error Detection in 
Computer Software", Proceedings of the American 

Ins titute of Indu strial Engi ne ers Natio nal Beeting, v. 

1, p. 1-6, May 21, 1975. 

4. Schneidewind, N. F. , "An Approach to Software 

Reliability Prediction and Quality Control", AF IP S 
Con fe ren ce P ro ceedings, v. 41, Part II, Fall Joint 
Computer Conference, p. 837-838, 1972. 

5. Baker, F. T., "Chief Programmer Team Management of 

Production Programming", IBM Sy st ems Jour nal, v. 11, 

No. 1, p. 65-66, 1972. 

6. Rizza, J. B. and Hacker, D., "Quality Assurance 

Inspection and Test Tools - An Application", 
Pro ce ed in gs of a Workshop on Currently Ava il able 

P r ogr am Testing Tools, v. 1, p. 9-10, April 1975. 

7. Howden, W. E., "Systems for Automating the Generation 
of Program Test Data", Proceedings of a Workshop on 
Currently Available Program Testing Tools, v. 1, p. 
37-39, April 1975. 


134 













Von Alven, W. H., (ed), Reliability Engineering, p. 
155-156, ARINC Research Corporation, Prentice Hall, 
1964. 

Stucki, L. G., "Automatic Generation of Self-Metric 
Software", Record of _197 3 IEEE Symposium on Computer 
Software Re li ability, v. New York City, p. 94-100, 
April 30 - May 2, 1973. 

McCormick, J. M„ and Salvadori, M. G., Num er ical 
Meth ods in F ORT RAN, p. 290-295, Prentice Hall, 1964. 










INITIAL DISTRIBUTION LIST 


4. 


♦ 


<• 

No. 


t 

% 


1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 

3. Department Chairman, Code 72 
Computer Science Group 
Naval Postgraduate School 
Monterey, California 93940 

4. Professor Norman F. Schneidewind, Code 32B 
Department of Operations Research 

and Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93940 

5. LT Thomas F. Green, USN 
3246 Fenelon Street 

San Diego, California 92106 


r 


r 


Copies 

2 

2 

1 

1 


136 




7 


r 


>4 JAM 76 '2 6 0(1(1 

22 OCT 76 2 3 6 16 * 


160963 

Thesis 

6742 Green 

c.2 Software error de¬ 

tection model. 

i 4 JAM 76 '*'24044 

22 OCT 76 2 3 6 16 l 


160963 

Thesis 

G742 Green 

c.2 Software error de¬ 

tection model. 





















